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ABSTRACT 



Measurement of the Transverse Single-Spin Asymmetry for 
Mid-rapidity Production of Neutral Pions in Polarized p+p 
Collisions at 200 GeV Center-of-Mass Energy 

Christine Angela Aidala 

The spin structure of the proton has revealed itself to be extremely complex 
and is an area of ongoing research. The Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory (BNL) inaugurated its operation as the first 
polarized-proton collider during the 2001-2002 run, marking the beginning of a new 
era in the study of proton spin structure. 

From the data collected in this run, the PHENIX experiment measured the 
transverse single-spin asymmetry (A N ) for neutral pion production at x F ^0.0 over 
a transverse momentum range of 1 to 5 GeV/c from polarized proton-proton inter- 
actions at a center-of-mass energy (y/s) of 200 GeV and found it to be zero within 
a few percent. Interest in these measurements arises from the observation of large 
(~30%) transverse single-spin asymmetries in p+p^— yjr+X at forward angles by the 
E704 collaboration at Fermilab (y/s = 19.4 GeV), found by the STAR and BRAHMS 
experiments to persist at RHIC energies, as well as single-spin, azimuthal asymme- 
tries observed recently in semi-inclusive deep-inelastic scattering experiments. Such 
large asymmetries were initially surprising because at leading order, perturbative 
quantum chromodynamics (pQCD) predicted only small effects. 

Several possible origins of these large asymmetries have been proposed. Despite 
great theoretical progress in recent years, no single, clear formalism has emerged 
in which to interpret the available data. Further theoretical work and a variety 
of additional experimental measurements will be necessary to understand current 
results and elucidate the transverse spin structure of the proton. 



For technical reasons, there have been minor textual modifications to the Table of 
Contents, List of Figures, and List of Tables in the present online version with 
respect to the final version of this thesis on record with Columbia University. 
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Chapter 1 
Introduction 

1.1 Proton structure 

The proton, together with the neutron and electron, is one of the basic building 
blocks of everyday atomic matter. Far from the point particle it was once believed 
to be, it has proven to be an extremely complex entity, and more than 80 years after 
it was discovered in the first decades of the twentieth century, the composition of 
the proton is still not completely understood. A very rich structure has gradually 
been uncovered over the past 40 years of research, with the appropriate description 
depending on the energy scale at which the proton is probed. The composition of 
the proton is now described by partons, including quarks and the particles carrying 
the force that binds them, gluons. More specifically, quarks can be categorized as 
either valence or sea. In the simplest composite model of the proton, it can be 
viewed as three valence quarks, each carrying 1/3 of the proton's linear momentum. 
In reality, measurements have demonstrated that there is a multitude of gluons and 
sea quark- ant iquark pairs present as well, each carrying generally a small fraction of 
the proton's momentum but with a large summed momentum contribution overall. 

1.2 The "proton spin crisis" 

In the naive model of the proton as simply three valence quarks of spin each, one 
might expect the proton's spin to be the straightforward sum of two parallel quark 
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spins and one antiparallel. However, in the late 1980's it was discovered that in fact 
only a small fraction of the proton's spin, less than 30%, was carried by quarks. This 
revelation, surprising at the time, came to be known as the "proton spin crisis." In 
retrospect, considering the complex linear- momentum structure of the proton, it is 
reasonable to expect a complex angular-momentum structure as well. Not only the 
spin of the partons is involved, but also their orbital angular momentum. 

Due to the fact that spatial rotations and Lorentz boosts do not commute, po- 
larized proton structure must be considered separately for a proton with spin vector 
parallel to or perpendicular to its (linear) momentum. This difference between 
the longitudinal and transverse spin structure adds further complexity to the prob- 
lem. As will be discussed, significant progress has been made in understanding 
the longitudinal spin structure of the proton, while transverse structure remains a 
largely open field. Despite impressive advances in just the past few years, numerous 
additional experimental measurements will need to go hand-in-hand with further 
theoretical investigation in order to elucidate the transverse spin structure of the 
proton. 

1.3 Studying proton spin structure at the Rela- 
tivistic Heavy Ion Collider 

Resolution of the proton spin crisis, in particular determination of contributions from 
sea quarks and gluons, remains the goal of extensive ongoing study. The PHENIX 
experiment at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory is in a unique position to make significant contributions to improve our 
understanding of the origin of the proton's spin. 

RHIC is the most versatile hadron collider in the world. It is capable of colliding 
heavy ions with energies as high as y/s = 200 GeV per colliding nucleon pair and 
polarized protons anywhere from 50 to 500 GeV, as well as different species in the 
two beams. In the first five years of running, RHIC has provided gold collisions at 
four different energies, copper collisions at three energies, deuteron-gold collisions, 
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and polarized-proton collisions. The flexibility of RHIC allows for a very diverse 
physics program. The heavy ion physics program investigates strongly-interacting 
matter at extreme temperatures and energy densities, seeking to create and study 
the properties of a state of matter known as the quark-gluon plasma (QGP). 

The polarized proton program seeks a better understanding of the proton's spin 
structure, in particular contributions to its longitudinal spin structure from the 
gluons and sea quarks. By studying the proton using hadronic collisions rather 
than electromagnetic probes, which do not couple to the electromagnetically neutral 
gluon, RHIC experiments may directly observe gluon-scattering processes. As a 
collider, RHIC can provide collisions at much higher energy than can be achieved in 
fixed-target measurements. As a result hard processes, describable by perturbative 
quantum chromodynamics (pQCD), can be studied, and new probes such as W 
bosons will eventually become available. 

1.4 Aims and outline of this thesis 

This thesis aims to motivate the study of the structure of the proton, in particular 
the polarized structure, as a fundamental question in QCD. It also seeks to describe 
how spin-dependent observables can be and have already been measured at RHIC, 
with focus on a transverse single-spin asymmetry (SSA) measurement, providing 
information on the transverse spin structure of the proton. 

This thesis will present a review of proton structure, both unpolarized and polar- 
ized, in terms of history and the current status. An overview of pQCD as applicable 
to proton-proton collisions at RHIC will be given. The RHIC polarized-proton ac- 
celerator complex and the PHENIX experiment and detector will be described. A 
measurement of the transverse single-spin asymmetry of neutral pions will be pre- 
sented and discussed. Finally, the prospects for future measurements to shed further 
light on the transverse spin structure of the proton will be explored. 
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Chapter 2 
Nucleon structure 

2.1 Unpolarized nucleon structure 

2.1.1 Elastic structure 

2.1.1.1 Magnetic moments 

The first evidence of proton substructure came from a measurement of its magnetic 
moment in 1933 by Esterman, Frisch, and Stern. It was found to be anomalously 
large and is now known to be approximately 2.79 times the Dirac magnetic moment, 
given by ~Ji p = -A- S , for a point-like spin-| particle of the same mass. The anoma- 
lous magnetic moment of the proton is now understood in terms of its valence quark 
structure (see Section 12 . 1 .21 below) and can be given by \x p = |(4/x u — fi d ), where \i u 
and fid are the magnetic moments of the up and down valence quarks, respectively. 

Similarly, the magnetic moment of the neutron was also found to be anomalous 
by Esterman and Stern in 1934 and is now understood in terms of its own valence 
quark structure. 

2.1.1.2 Form factors 

Charge and current distributions within the nucleon can be described by electro- 
magnetic form factors, measurable via elastic electron-proton scattering. Viewed in 
a particular frame known as the Breit frame (p final — ~~ P initial f° r the proton), the 
form factors Ge and Gm are proportional to the Fourier transforms of the charge 
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and magnetization distributions, respectively. 

The cross section for elastic electron-proton scattering can be expressed in terms 
of the form factors, as given in Eq. 12. 1[ 



where r = —q 2 /4M 2 , q is the four- momentum transfer in the scattering, M is the 
proton mass, a is the fine structure constant, 9 is the electron scattering angle in the 
laboratory frame, and E and E' are the incident and scattered electron energies. 
The electric and magnetic form factors can be determined either via differential 
cross section measurements of unpolarized e + p scattering or via measurement of 
the recoil proton polarization in the reaction e + p — > e + p , where the arrows 
indicate polarization. For a relatively recent summary of proton electromagnetic 
form factor measurements performed at Jefferson Lab, see |129j . 

2.1.2 Inelastic structure: QCD and the quark-parton model 

In the 1960's deep-inelastic lepton-nucleon scattering (DIS) experiments at SLAC, 
analogous to the famous Rutherford scattering experiment that led to the discovery 
of the atom's hard core, found that protons also had "hard" subcomponents |62|I73]. 
These hard subcomponents came to be known as partons. It took some time before 
the experimentally observed partons inside the proton were identified as the so-called 
"quarks," which had been theoretically hypothesized based on hadron spectroscopy 
data as part of the "Eightfold Way" by Gell-Mann and Ne'eman independently in 
the early 1960's (HHI- But eventually the quark-parton model of the proton came 
into being. As experimental work progressed and higher-energy lepton beams were 
used as probes, the proton came to reveal a much more intricate structure than that 
of the three so-called "valence" quarks. These other subcomponents are now known 
to be sea quarks and gluons. 

The experimental and theoretical work in the 1960's and 1970's regarding hadronic 
interactions and structure led to the development of the theory of quantum chro- 
modynamics (QCD), describing the behavior of the strong force. A central concept 
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of QCD is that of asymptotic freedom. While quarks are strongly bound at distance 
scales larger than a typical hadron radius (r « 1CT 15 m), at shorter distances they 
behave as nearly free. 

2.1.2.1 Structure functions 

The nucleon structure functions describe the inelastic structure of the proton and 
neutron, probed principally via DIS. The double differential cross section for inelastic 
cross section for electron-proton scattering can be expressed as in Eq. 12.2^ 



d 2 o . Aa 2 E' 2 [~ . 2 , 2 9 nTr/ , 2 , . 2 9 
llab = — W 2 (v, q 2 ) cos 2 - + 2W x (v, q 2 ) sin 2 - 



(2.2) 



dE'dfl' q A 

in which q, a, 9, and E' are as in Eq. 12.11 v — (p ■ q)/M with p being the initial 
nucleon four-momentum, and W\ and W2 are the proton structure functions. Note 
the similarity of Eq. 12.21 to Eq. 12.11 the cross section for elastic electron-proton 
scattering, with the structure functions playing the role of the form factors. 

It is common to express the proton structure functions slightly differently, as 
given in Eq. 12.31 

Fx = MWt 

F 2 = uW 2 (2.3) 

Fx and F 2 can be written as functions of Q 2 = —q 2 > and the dimensionless 
variable x = Q 2 /2p ■ q. They can be related to the cross sections for scattering of 
transversely and longitudinally polarized virtual photons off of the proton, ot and 
<7l, as given in Eq. 12.41 

o L oc (P--F1) (2.4) 



2x 

The total virtual photon-proton cross section is the sum of these components and 
is proportional only to F 2 . 

In 1969 Bjorken predicted that at large Q 2 , scattering off of "point-like" subcom- 
ponents that were approximately free in the proton would lead to proton structure 
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functions with no Q 2 dependence for a given value of x, i.e. that scale with just this 
single, dimensionless variable jHJ. That is, assuming point-like constituents of the 
proton, at large Q 2 , inelastic electron- proton scattering could be viewed as elastic 
scattering of an electron off of a hard, point-like particle within the proton. 

The experiments performed at SLAC mentioned above jH2 EH] discovered the 
scaling behavior predicted by Bjorken. The structure functions they measured had 
very little explicit dependence on Q 2 and could in fact be written simply as functions 
of x. 

In Figure 12.11 showing world DIS data for F 2 of the proton, one sees that for 
x > 0.02, F2 is nearly flat in Q 2 . This indicates that the hard subcomponents being 
probed are approximately free. The early measurements included here provided 
evidence leading to development of the concept of asymptotic freedom in QCD. 
Note, however, the scaling violations observed at low x and low Q 2 . These are now 
understood in terms of gluon radiation emitted by the parton prior to the hard 
scattering. See Section 12.1.2.21 for discussion of the relation of scaling violations to 
the gluon distribution function. 

In 1969 Callan and Gross predicted that for spin-| charged components within 
the nucleon, the scaling structure functions would be related as given in Eq. 12.51 
known as the Callan-Gross relation [77] . 

F 2 (x) = 2xF 1 (x) (2.5) 

Experimental confirmation of the Callan-Gross relation came from SLAC in the late 
1970's [HI] and thus provided strong evidence for the spin-| nature of what are now 
known to be quarks. 

2.1.2.2 Parton distribution functions 

Feynman's introduction of the quark-parton model (QPM) in 1969 jHH| offered a 
relatively intuitive explanation of Bjorken scaling. The virtual photon in DIS could 
be viewed as scattering elastically off of a collection of hard partons within the 
proton; the DIS cross section is the incoherent sum of the individual cross sections. 
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Figure 2.1: World DIS data on the unpolarized structure function of the proton, 
shown compared to next-to-leading order (NLO) pQCD fits. 
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As the proton momentum approaches infinity, x can be seen as the fraction of the 
proton's linear momentum carried by the parton. The size of the cross section for 
scattering off of a particular parton is proportional to the probability, q(x), of hitting 
a quark of flavor q carrying momentum fraction x of the proton. q(x) is known as 
a parton distribution function (pdf). The scaling structure functions, F\(x) and 
F 2 (x), can then be viewed as representing the probability of scattering off of a 
parton within the proton carrying momentum fraction x. They can be expressed 
in terms of the pdf's for different quark flavors as follows, where indicates the 
electromagnetic charge of the quark of flavor i. 



Including all partons in the proton and not only the charged particles, which 
couple electromagnetically and can be probed directly by DIS, a momentum sum 
rule is obtained (Eq. 12. 7|) . 



It has been found experimentally that approximately half of the proton's total mo- 
mentum is carried by (electromagnetically neutral) gluons, which are present as 
mediators of the strong interactions among the quarks in the nucleon and dominate 
at low x values. In the Bjorken-scaling regime, at large x, the interactions among 
the quarks, i.e. the gluons, do not play a significant role. While DIS experiments 
cannot directly probe the gluon distribution function, g, the scaling violations and 
the evolution of the structure functions in Q 2 provide information on the gluon. 
The gluon distribution can be obtained from the logarithmic scaling violations of 
the structure function F 2 , as given in Eq. 12.81 







(2.7) 



g oc dF 2 /d{\nQ 2 ) 



(2.8) 
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See Section 12.1.41 for a discussion of the current status of pdf measurements. 

2.1.3 pQCD, factorization, and universality 

Performing calculations in QCD presents a number of challenges that quantum 
e/ectrodynamic (QED) calculations do not. In QCD the force carriers themselves 
are charged; gluons carry color charge, whereas photons, the force transmitters in 
QED, are electrically neutral. Contributions from higher-order Feynman scattering 
diagrams in QED, i.e. higher powers in the electromagnetic coupling constant, a ~ 
1/137, representing additional lepton-photon vertices, quickly become negligibly 
small, due to the fact that a is much less than one. An analogous expansion in 
QCD is only possible in the regime where the coupling, a s (Q 2 ), is small, which is 
generally the case for processes involving a large momentum transfer. Perturbative 
QCD (pQCD) is the calculation technique used in this kinematic regime. 

A hadron can be viewed as a collection of free, massless partons with parallel 
momenta. The collinear factorization theorem in pQCD starts from this assumption 
of collinearity of the partons and hadrons, i.e. no transverse momentum of the par- 
tons in the proton with respect to the initial proton momentum, and no transverse 
momentum of the final-state hadron with respect to the scattered parton momen- 
tum. The collinear factorization theorem separates cross sections for hard-scattering 
processes into parts that are soft, or non-perturbative, and hard, or perturbative, 
in a self-consistent way. The soft components, pdf's and fragmentation functions 
(FF's), must be obtained from experimental measurements. The hard components, 
partonic hard-scattering cross sections, are directly calculable in pQCD. Parton dis- 
tribution functions, discussed above in Section 12.1. 2. 2\ correspond to the probability 
of striking a particular parton carrying momentum fraction x of the proton; FF's 
represent the probability of the scattered parton fragmenting into a particular final- 
state hadron, as a function of the fraction z of the scattered parton's momentum 
passed along to the final-state hadron. 

The factorization theorem was developed and proven over the course of the late 
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1970's to the mid-1980's. Early work can be found in [HE1 H3 EH EQ|; 
complete proofs are available in IHH} 11241 . The factorized cross section for 
hard scattering in hadron-hadron collisions (A + B — > C) is given by Eq. 12. 9| 

rfa = ^ f a (x a ,n f ) <g> / 6 (z 6 , ///) <g> da ab (x a , x b , z c ,n f ,n r ) <g> Z)f (z c , ////) (2.9) 

abc 

in which / a (/&) is the density of parton a (b) in hadron A (B), is the fragmenta- 
tion function of parton c into hadron C, and /if and //// are arbitrary scales known 
as factorization scales, which can be thought of as the amount of parton radiation 
incorporated into the pdf 's and FF, or as the separation scale chosen to distinguish 
between the hard and soft components of the cross section. While the scales chosen 
are arbitrary, they must be chosen consistently between the soft components and the 
partonic hard-scattering cross section, da a b. The partonic cross section depends on 
an additional arbitrary scale, the renormalization scale, which controls the running 
of the strong coupling, a s . It is common practice to set the factorization and renor- 
malization scales to be equal; a typical value chosen is one close to the momentum 
transfer of the process. 

The principle of universality in conjunction with the factorization theorem makes 
the formalism of pQCD extremely powerful. The principle of universality states that 
pdf's and FF's are the same regardless of the scattering processes involved. Univer- 
sality implies the dominance at high momentum transfer of leading- twist (twist-two) 
contributions, with interactions only between the two hard-scattering partons. A 
higher-twist calculation takes into account the exchange of additional gluons be- 
tween the hard-scattering partons and the nucleon remnants. The twist expansion 
is in successive powers of 1/Q 2 ; therefore, higher-twist contributions are suppressed 
for processes with large momentum transfer. Because of universality, pdf's and FF's 
can be measured in the environment which allows the most accurate determination 
and then utilized as input for pQCD calculations in other processes. For example 
quark distribution functions can be measured in DIS experiments, in which the kine- 
matics and thus the probed x values are straightforward to understand, then utilized 
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in calculations for the more complicated environment of hadron-hadron collisions. 
FF's are most easily measured in e + + e~ collisions because the four-momenta of the 
outgoing quarks are well known. Decades of comparison between experimental cross 
section measurements and pQCD have provided a testing ground for the assumption 
of universality, and by now it is a well established and accepted principle. 

While pdf's are not calculable in pQCD and are typically obtained from exper- 
iment, they are in principle calculable using other theoretical techniques such as 
lattice QCD. For recent calculations of structure and distribution functions on the 
lattice, sec [128J and references therein. 

Over the course of the 1970's, a formalism emerged in which it was possible 
to take a measurement of a pdf at a particular value of x and Q 2 and predict 
the pdf at the same x but different Q 2 . This formalism is known as DGLAP, 
acknowledging important contributions from Dokshitzer, Gribov, Lipatov, Altarelli, 
and Parisi |101| I119[ l98| 130] . DGLAP has been essential to the relevant application 
of factorized pQCD, which requires as input experimentally measured pdf's, which 
necessarily are available at only a finite set of x and Q 2 values. Using DGLAP, 
calculations can be done for any Q 2 value desired. 

There is additionally a prescription for evolution of measured pdf's to different 
values of x for fixed Q 2 , formulated by Balitsky, Fadin, Kuraev, and Lipatov, also 
in the 1970's [EH I116| I117j . The BFKL technique has demonstrated itself to be 
similarly useful in performing calculations at desired values of x and Q 2 . 

Factorized pQCD has proven to be a valuable and successful theoretical technique 
for many years now. Its applicability to measurements in p + p collisions at RHIC 
will be discussed in Chapter El 

2.1.4 Current status of the unpolarized structure of the pro- 
ton 

A long history of experiments, in particular DIS experiments at SLAC, CERN, and 
DESY, has measured the unpolarized structure of the proton well. Gluons have 
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been found to play an important role, carrying approximately 50% of the proton's 
momentum. A comprehensive review of the contributions to unpolarized proton 
structure made by the experiments at the HERA electron-proton collider at DESY 
is given in A thorough review of nucleon structure functions and pdf 's is available 
in [EH]. 

The wealth of accumulated data regarding the unpolarized structure of the pro- 
ton is evident in Figure 12.11 Measured x values range from deep into the sea at 
x = 3.2 x 1CT 5 to well into the valence region at x — 0.65, with Q 2 values as high 
as 10 4 GeV 2 . Periodic efforts have been made to examine all data available and 
perform a global analysis in order to obtain the best-fit pdf's; see for example |131j . 

In the 1990's it was discovered by the NMC experiment at CERN that there is a 
flavor asymmetry in the unpolarized sea of light quarks in the proton [SSI US] • There 
is a significant excess of d with respect to u. Although no known symmetry requires 
d/u = 1, the experimental result was unexpected. It has since been confirmed by 
other experiments at CERN [HH], Fermilab fTTT^irTIK] . and DESY jij but is still not 
well understood. 

2.2 Polarized proton structure 
2.2.1 Historical overview 

For many years it was assumed that the proton's spin of \fi was due to the spins of 
the three spin-| valence quarks, with two oriented in one direction and one in the 
other. In the late 1980's, however, the EMC experiment at CERN |46| l4Tj discovered 
that only approximately 13 ± 16% of the proton's spin was due to the spin of the 
quarks. This surprising result became known as the "proton spin crisis." With so 
little of the proton's spin coming from the total quark spin (AE), the remainder 
is expected to come from gluon spin contributions (Ag) and the orbital angular 
momentum (OAM) of both quarks and gluons (L g+q ), as indicated in the spin sum 
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rule given in Eq. 12.101 which is valid in the infinite momentum frame. 

~ = ^AS + As + L g+g (2.10) 

Experimental work following the EMC discovery, mostly exploiting DIS, has con- 
tinued to explore this problem for more than 15 years, yet there remains much to 
be understood. In particular, the magnitude and even sign of the gluon spin contri- 
bution to the spin of the proton remains to be determined, the flavor breakdown of 
the sea quark spin contributions is largely unknown, and no definitive experimental 
technique with which to access OAM directly has yet been proposed. 

2.2.2 Polarized structure functions 

Similar to the unpolarized structure functions discussed in Section 12.1.2.11 spin- 
dependent structure functions can also be defined. The difference in cross sections 
for deep-inelastic scattering of leptons polarized antiparallel and parallel to the spin 
of the target proton can be written as in Eq. 12.111 



d 2 (T + ~ d 2 a ++ Ann 2 

1QH V -WT V - ^\M(E + E>. M ^)-Q^m (2.1!) 

in which the kinematic variables are defined as in Equations 12.11 and 12.21 and G\ 
and G2 represent polarized structure functions of the proton. In the Bjorken scaling 
limit of large Q 2 and u, these structure functions depend only on x and can be given 
as in Eq. EH 



g x {x) = M 2 uG 1 (u,Q 2 ) 

g 2 (x) = Mv 2 G 2 {y,Q 2 ) (2.12) 

gi(x) can be viewed as the difference in probability of scattering off of a parton 
carrying momentum fraction x of the proton with parton helicity antiparallel versus 
parallel to the proton spin. Figure EZ3 shows the world DIS data for g\ as a function 
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Figure 2.2: World DIS data on the polarized structure function of the proton. 



of Q 2 ; the scaling behavior can be seen over most of the kinematic range that has 
been measured, which is limited compared to the kinematic range over which F2 has 
been measured (refer back to Figure 12. lj) . 
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2.2.3 Polarized parton distribution functions 

Polarized pdf's, or helicity distribution functions, denoted Af, represent the differ- 
ence in probability of scattering off of a parton / with its spin vector parallel versus 
antiparallel to the proton's spin, in the case of longitudinal polarization, i.e. po- 
larization along the direction of proton motion. The transverse spin structure of 
the proton is discussed separately in Section 12.2.51 Helicity distributions can be 
obtained for example from global fits to DIS measurements of g±(x) for the proton 
and neutron. It is well accepted that polarized pdf's can be used as input to calcu- 
lations in pQCD in a similar fashion to unpolarized pdf's. Arguments for collinear 
factorization involving spin- dependent processes are given in [75] . 

2.2.4 Current status of the longitudinally pol. structure of 
the proton 

There are multiple ongoing experiments making measurements to study the po- 
larized structure of the nucleon. The HERMES experiment at DESY makes use 
of the longitudinally polarized electron (or positron) beam at the HERA collider 
and performs DIS measurements on nucleon targets which can be longitudinally or 
transversely polarized. The COMPASS experiment at CERN also performs spin 
measurements through DIS, using a polarized muon beam from pion decays on a 
polarized fixed target. The PHENIX and STAR experiments at RHIC study po- 
larized proton-proton collisions and have the ability to choose either longitudinal 
or transverse beam polarization, while the BRAHMS experiment at RHIC can only 
study transverse spin physics (see Section l4~2|) . 

In 2004 HERMES published its final longitudinal results, the first five-flavor fit 
for quark polarizations measured in DIS [21]. Their results yield Au/u > and 
increasing with x, Ad/d < 0, and polarizations consistent with zero for u, d, and s 
quarks. In contrast to the unpolarized case, the polarized light-quark sea does not 
appear to have a significant flavor asymmetry. 

Various different groups have performed global fits to the world data available 
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from DIS experiments on the helicity structure of the proton jHSl HOOj I104j . Fig- 
ure 12.31 depicts the results of several such fits for valence quarks; Figure 12.41 for 
gluons and sea quarks. Both figures are taken from |l()4j : see the reference for de- 
tails regarding the different curves. The best constraints are naturally available for 
valence quarks. As can be seen, the polarization for valence up quarks is significant 
and positive, while for valence down quarks it is significant and negative, consistent 
with the HERMES results mentioned above, which were not entirely available at 
the time of these fits. A global analysis of the sea quarks yields a negative polariza- 
tion of smaller magnitude than for valence down quarks and with a larger relative 
uncertainty. 

As mentioned above for the unpolarized case, DIS experiments cannot access the 
gluon directly because it does not couple electromagnetically. Consequently, Ag can 
only be inferred through scaling violations of the g\ polarized structure function, as 
given in Eq. 12.131 

Ag oc dg l /d(\nQ 2 ) (2.13) 

or via di-hadron or heavy flavor production, both of which provide some sensitivity 
to the gluon. The magnitude of Ag remains almost completely unknown, and its 
sign is not yet clear, as can be seen in Figure 12.41 A recent comprehensive review 
of the longitudinal spin structure of the proton can be found in [5^] . 

2.2.4.1 Recent results from RHIC 

From the 2003 polarized proton run at RHIC, a measurement of the longitudinal 
double-spin asymmetry (All) of neutral pions at mid-rapidity has been made by 
PHENIX for a transverse-momentum (p?) range of 1 to 5 GeV/c. Pion pro- 
duction in this kinematic region is due mainly to gluon-gluon and gluon-quark scat- 
tering; thus, this measurement is sensitive to the polarized gluon distribution. (See 
Figure 17.11 in Chapter [7| indicating the relative contributions of different partonic 
scattering processes to tt production as a function of px-) The asymmetry is shown 
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Figure 2.3: Global fits to world data for the polarized pdf 's of valence quarks, taken 
from jlU4j . 

in Figure 12.51 The theoretical curves in the figure represent next-to-leading-order 
(NLO) pQCD calculations using two different assumptions for Ag |1(JU[ 111)7] . GRSV- 
std takes Ag as the value that best fits the world DIS data, and GRSV-max takes 
Ag to be equal to the unpolarized gluon distribution at a scale of Q 2 = 0.6 GeV 2 . 

Due to the significant contribution of gluon-gluon scattering and the isospin 
symmetry of the neutral pion, its double-spin asymmetry is largely insensitive to 
the sign of the polarized gluon distribution. Future results for All of positive and 
negative pions for px > 5 GeV/c will provide independent measurements of Ag 
and allow determination of its sign. See Figures 12.61 and 12.71 for expected charged 
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Figure 2.4: Global fits to world data for the polarized pdf 's of gluons and sea quarks, 
taken from 



pion asymmetries as a function of pt for different-sign polarized gluon distributions, 
calculated by M. Stratmann. PHENIX expects to be able to make a significant mea- 
surement of charged pions in the next long polarized proton run at RHIC, expected 
to occur in 2006 or 2007. All of direct photon production, measurable on a slightly 
longer time scale, will provide a clean measurement of both the magnitude and sign 
of Ag, as discussed in Section POl 
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Figure 2.5: Longitudinal double-spin asymmetry for neutral pion production at 
PHENIX, compared to predicted asymmetries assuming various values of Ag (see 
text), taken from [T7] . 

2.2.5 Understanding the transverse spin structure of the 
proton 

2.2.5.1 Longitudinal versus transverse spin structure of the proton 

The transverse spin structure of the proton cannot be determined from its longitu- 
dinal spin structure. A simple explanation for this fact is the non-commutation of 
Lorentz boosts and spatial rotations. Differences between the transverse and lon- 
gitudinal polarization structure of the nucleon provide insight into the relativistic 
nature of partons bound within the nucleon. A concise description of transverse spin 
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Figure 2.6: Predicted All of positive pions for various Ag assumptions, from M. 
Stratmann. 

0.08 
0.06 
0.04 
0.02 


-0.02 
-0.04 

2 4 6 8 10 12 14 

p T (GeV/c) 



Figure 2.7: Predicted All of positive pions for various Ag assumptions, from M. 
Stratmann. 

structure functions and their relation to the longitudinal structure of the proton can 
be found in 106J. 

Inelastic and elastic scattering can be related by the optical theorem. Thus 
inelastic proton-proton or quark-proton scattering can be considered in terms of 
elastic quark-proton scattering. In the elastic scattering of two spin-| particles, 
there are three different possibilities for the initial- and final-state helicities. The 
particles may start and end with the same helicity (++ — > ++), start and end 
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with opposite helicities (H > H — ), or start with opposite helicities and change 

helicities in the scattering (H > — h). Linear combinations of these three helicity 

configurations in the scattering can be formed, corresponding to the momentum (q), 
helicity (Ag), and transversity (Sq) distribution functions given in Eq. 12.141 

q : (++ -> ++) + (+- ->+-) 
Ag : (++ -> ++) - (+ ► +-) 

<5g : (+--►-+) (2.14) 

Transversity is therefore a chiral-odd, or "helicity- flip," distribution. In a trans- 
verse basis, it represents the difference in probability of scattering off a quark with 
transverse spin parallel versus antiparallel to a transversely polarized proton. This 
interpretation is directly analogous to the meaning of helicity distributions in a he- 
licity basis. The transversity distribution was first discussed in The Soffer 
bound, given by Eq. I2.15[ relates the transversity, helicity, and momentum distri- 
butions of the nucleon |142j . 

\25q(x)\ < q(x) + Aq(x) (2.15) 

Analogous to the longitudinal spin sum rule (Eq. 12. 10]) . there is also a transverse 
spin sum rule [54J, given in Eq. 12.161 in which L$ T is the transverse component of 
the partonic orbital angular momentum. 

\ = \ E / dxS( la(x,Q 2 ) + Yl ( L S T )a(Q 2 ) (2-16) 
a=g,q a=q,q,g 

Note the absence of a gluon spin contribution; there is no transversity distribution 
for gluons at leading twist because there is no mechanism to flip the helicity of 
(spin-1) gluons in the scattering. 

While transversity is a non-perturbative object as any other pdf, the nucleon 
tensor charge, related to transversity (J^ dx(8q(x) — 5q(x))), can be calculated in 
lattice QCD. For a description of recent work, see |128j . 
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In pQCD, chiral-odd functions must appear in pairs because hard scattering 
processes conserve helicity. One possibility is to look for observables that represent 
the convolution of two transversity distributions, i.e. double transverse- spin asym- 
metries, Att- Another possibility is to convolute transversity with a chiral-odd 
fragmentation function, one example of which will be discussed below. 

2.2.5.2 Observation of large transverse single-spin asymmetries 

The measurement of transverse single-spin asymmetries (SSA's), for example in 
proton-proton collisions or DIS, represents one way of probing the quark and gluon 
structure of transversely polarized nucleons and is the approach exploited for the 
measurement in this thesis. Interest in these measurements is heightened by the 
large transverse SSA's observed in spin-dependent proton-proton scattering exper- 
iments spanning a wide range of energies. The experimental observation of large 
asymmetries, with the first measurements being in the late 1970's, was initially a 
surprise. The leading-twist pQCD expectation was that transverse SSA's should be 
suppressed as ^^ £ , where m q is the quark mass |112j . 

Striking asymmetries were seen at a number of spin-dependent p + p scattering 
experiments at energies ranging from y/s = 5 — 10 GeV. Asymmetries approaching 
30% were observed in inclusive pion production at large Feynman-x (xp = 2pp/\/s, 
where pi is the component of particle momentum in the beam direction) and pp 
up to 2 GeV/c |94[ l2*9*j . In a different kinematic region at mid-rapidity and large 
xt = 2px/y/s, asymmetries were also observed in inclusive 7r° and ir + production but 
not in 7r~ production |47 H I42j. At a higher center-of-mass energy of 19.4 GeV 
where pQCD may be applicable, asymmetries at large xp persisted [310]; however, 
the asymmetry in 7r° production at mid-rapidity at this energy was found to be 
zero up to pt of 4 GeV/c jH]. Non-zero transverse single-spin asymmetries were also 
observed in semi-inclusive DIS jTHl EI3 H2j • 

The results for neutral and charged pions at high xp from the E704 experiment 
at Fermilab [B] are shown in Figure 12.81 The observed asymmetries are strikingly 
large, reaching a magnitude of ~ 40% for charged pions at xp ~ 0.8. There is a clear 
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Figure 2.8: Transverse SSA of high-x^ neutral and charged pions at ^/s = 19.4 GeV, 
taken from j^j. 

sign dependence of the asymmetry on the pion charge, with A]^ > and A* N < 
and both of approximately equal magnitude and exhibiting the same dependence 
onij?. The tt° asymmetry is also positive but with smaller magnitude. 

Generally, transverse SSA's can be described by spin-momentum correlations 
having the form S ■ (pi x p%), in which there are different possibilities for the 
spin and momentum three-vectors. Three different mechanisms, originating from 
different spin-momentum correlations, have been studied extensively as the possible 
origin of transverse SSA's in high-energy hadron collisions. 



1. Transversity distributions can give rise to SSA's in combination with spin- 
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dependent, chiral-odd FF's, which serve as analyzers for the transverse spin 
of the struck quark. The Collins function jTS] is an example of such a FF. 

2. Quark and gluon distributions that are asymmetric in the intrinsic transverse 
momentum of the parton within the proton, k?, can lead to SSA's. This idea 
was first suggested by Sivers |14()j . The Sivers pdf can exist both for quarks 
and gluons, and a relation to orbital angular momentum of partons in the 
nucleon has been suggested |14U| ITH] . 

3. Interference between quark and gluon fields in the initial or final state can 
generate SSA's [1321 ITTT] . 

2.2.5.3 Non-collinear pdf's and FF's 

The polarized pdf's and FF's relevant to collinear factorization given in Section l2.2.3l 
are integrated over all possible values of k^. In a pdf, kx represents the transverse 
momentum of the parton within the proton; in a FF, it indicates the transverse 
momentum of the fragmenting hadron with respect to the scattered quark, or jet 
axis. Transverse-momentum-dependent (TMD) pdf's relate naturally to the orbital 
angular momentum of partons within a proton; however, a precise understanding of 
this relation remains unclear. If one does not assume collinearity but rather takes 
/cy-dependent pdf's and FF's, a rich variety of new possibilities emerges. 

There are a total of eight leading-twist (twist-two) /c^-integrated and /^-dependent 
quark distribution functions, as shown in Figure 12.91 Only three of these pdf's are 
independent of kx- Parton distributions denoted by / indicate unpolarized quarks, g 
longitudinally polarized quarks, and h transversely polarized quarks. The subscript 
1 signifies leading twist; the subscripts L and T denote longitudinal and transverse 
proton polarization, respectively. The superscript _L indicates explicit dependence 
on transverse momenta with a non-contracted index, as described in 1.25.. Note 
that hyr is the transversity distribution, an alternative notation for Sq. The field 
of transverse spin physics has been plagued for many years by confusing and incon- 
sistent notation in the literature for relevant structure functions, pdf's, and FF's, 
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Figure 2.9: The eight leading-twist quark distribution functions, including collinear 
as well as /^-dependent distributions. 



as well as sign conventions for azimuthal angles in semi-inclusive DIS. One of the 
outcomes of the Transversity 2004 workshop in Trento, Italy was an examination 
and comparison of what existed in the literature and a set of recommended notation 
and sign conventions |49j . 

It should be pointed out that the factorization theorem has not been proven 
generally for the case of non-collinear partons. It has so far only been proven for 
the Drell-Yan process (q + q — ► £ + + £~) [HZ], with notable work also in e + + e~ — > 
2h + X and semi-inclusive DIS |831 185]. /cr-dependent factorization is therefore 
strictly speaking an assumption, albeit a well-accepted one in the field. Efforts are 
ongoing to establish the theoretical basis more firmly. Recent work considering gauge 
invariance in the cases of semi-inclusive DIS and Drell-Yan appears in |l()9[ll()8t l82j. 
Discussions of the universality and the evolution of TMD pdf's can be found in 
[HI EE22 E3 and [HEl 1103 HIS IMl, respectively. In general, the role of k T in 
hard-scattering processes is a vibrant area of research in pQCD. 



28 



CHAPTER 2. NUCLEON STRUCTURE 



2.2.5.4 Sivers effect 



One of the pdf's that emerges if is left unintegrated is known as the Sivers 
pdf, denoted as f^ T and first proposed by Sivers in 1989 |14()[ 1141] , It represents 
the spin-dependent asymmetry in the intrinsic &t of the (unpolarized) partons in a 
transversely polarized proton. The Sivers function was for some time believed to be 
forbidden because it is time-reversal odd (T-odd), but it was finally realized in 2002 
that final-state interactions via a soft gluon can create the necessary interference 
of amplitudes for the Sivers function to exist [71]. The T-odd nature of the Sivers 
function is now commonly referred to as "naive T-odd" to express the fact that it 
is not in fact forbidden in QCD. 

The Sivers function plays a central role in the phenomeno logical Sivers effect, 
which has its origin in correlations of the form S ■ (P x fcy), where S is the proton 
spin, P the proton momentum, and kr the intrinsic transverse parton momentum 
in the proton. In a simplistic picture of the Sivers effect, the transverse polarization 
of the proton can be viewed as originating from the orbital angular momentum of 
the partons. A spin-dependent final-state azimuthal asymmetry is then generated 
by preferential scattering off of the orbiting partons in the "front" or "back" of 
the proton, with scattering off of the front of a proton with spin up generating 
particle production preferentially to the left of the polarized beam. The mechanism 
responsible for scattering off of a particular "side" of the proton is not entirely 
clear. Additional study of the Sivers effect has been performed by Burkardt as well 
as Hwang [23 EE] • 

The absolute value of the Sivers function would provide a lower bound on parton 
OAM. The Sivers function was originally investigated for quarks; that for gluons was 
first discussed only in 2003 |138j . Asymmetry calculations based on the Sivers effect 
to describe the E704 data and other results can be found for example in [3E1IHH]- 
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2.2.5.5 Collins effect 

As mentioned above, chiral-odd functions in pQCD must come in pairs because 
helicity is conserved in hard scattering processes. Therefore transversity, as a chiral- 
odd distribution, needs to be convoluted with another chiral-odd function in order 
to be relevant in physical processes. For production of final-state hadrons, a chiral- 
odd fragmentation function is one possibility. Such a FF was proposed in the early 
1990's by Collins, Heppelmann, and Ladinsky [7EJ HE EE] • 

The Collins FF, denoted H^~, represents the correlation between the transverse 
polarization of the fragmenting quark and the orientation of the hadron production 
plane, given by S • ( k x Ph), in which S is the transverse polarization of the 
scattered quark, k is its three-momentum, and P h is the three- momentum of the 
final-state hadron. 

The Collins mechanism, a phenomenological mechanism incorporating the Collins 
FF, describes a spin-dependent azimuthal asymmetry in the distribution of hadrons 
within a jet. A relatively intuitive model of the Collins mechanism for the produc- 
tion of pseudoscalar mesons has been proposed by Artru et al. |45j . In this model, a 
transversely polarized quark is scattered out of a transversely polarized proton, with 
the probability of the direction of the scattered quark spin given by the transversity 
distribution. In order to produce a pseudoscalar (spin-0) meson such as a pion, the 
fragmenting quark must acquire an oppositely polarized (anti-) quark from the vac- 
uum. If the quark-antiquark pair from the vacuum is assumed to have a total spin 
angular momentum of 1, conservation of angular momentum requires one unit of 
orbital angular momentum in the opposite direction. This orbital angular momen- 
tum of the (anti-)quark from the vacuum which subsequently binds to the scattered 
quark then produces a preference in azimuthal direction in the production of the 
final-state pion. 

While it was originally believed that it was possible to explain the large transverse 
SSA's observed entirely in terms of the Collins effect (see for example 01]), recent 
work has suggested that this might not be the case |3Blll20j. As yet, no consensus 
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has been reached. 

2.2.5.6 Higher-twist effects 

It has been shown that higher-twist contributions and non-zero kx can produce the 
same effects in hard-scattering processes |130j . As such, there have been studies 
of how twist-three effects rather than TMD distributions can give rise to the large 
SSA's observed. 

Qiu and Sterman have examined higher-twist asymmetry contributions in the 
initial state system, i.e. interference between quark and gluon fields in the polarized 
proton [132J. Similar studies have been performed by Kanazawa and Koike for 
quark-gluon interference in the final state, a parton fragmenting to a hadron 
In the initial-state case, both chiral-even and chiral-odd components are possible. 
It is believed that a relation of the chiral-odd twist-three initial-state effect of [132 
to the Sivers mechanism may exist, but this is not completely understood. 

2.2.6 Recent experimental results in transverse spin physics 

RHIC data for polarized-proton collisions at ^/s = 200 GeV first became available in 
late 2001. In addition to the results from the PHENIX experiment presented in this 
thesis, transverse spin measurements have been made by the STAR and BRAHMS 
experiments at RHIC. At STAR it was discovered that large transverse SSA's persist 
even at RHIC energies, an order of magnitude higher than the energy for previous 
results jH]. They found asymmetries of up to ~ 30% in the forward production 
of neutral pions, as can be seen in Figure 12.101 The theoretical curves represent 
different fits to the E704 results discussed above 0, scaled in energy from 19.4 to 
200 GeV. The general agreement of the scaled fits with the 200-GeV data suggests 
that the asymmetries are generated by similar mechanisms at the two energies; 
however, further study is needed before any definitive conclusion can be made. 
STAR also has preliminary asymmetry results for the production of neutral pions in 
the backward direction with respect to the polarized beam, for —0.6 < xp < —0.2, 
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Figure 2.10: Transverse SSA of high-x^ neutral pions at ^/s = 200 GeV, taken from 
0. See text for more details. 

which are consistent with zero [123 . 

BRAHMS has preliminary results for the transverse SSA's of charged pions as 
well as protons |146j . The charged pion asymmetries are shown in Figure 12.111 
The charge dependence of the sign of the asymmetry clearly follows that observed 
by E704 [6J; however, the results are for low transverse-momentum values (p? < 
3 GeV/c), so pQCD-based interpretations such as the Sivers and Collins effects may 
not be applicable. The results are in reasonable agreement with extrapolations of 
initial-state twist-three calculations by Qiu and Sterman. 

The HERMES experiment at DESY has measured non-zero Collins and Sivers 
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Figure 2.11: Transverse SSA's of moderate-277 charged pions at -^s = 200 GeV, 
taken from |146j . The theoretical curve is an extrapolation of initial-state twist- 
three calculations by Qiu and Sterman. 

moments via semi-inclusive deep-inelastic scattering of positrons off of a hydrogen 
target [22] • The COMPASS experiment at CERN made a similar measurement 
with a muon beam on a deuteron target and found results consistent with zero 
[2*7] . The COMPASS observation is now understood to be because of asymmetry 
cancelations due to the isospin symmetry of the deuteron target at COMPASS. 
For an interpretation of the recent HERMES and COMPASS results and related 
predictions for transverse SSA's at RHIC, see [147]. 

The Collins FF for pions was recently measured by the BELLE e + + e~ annihila- 
tion experiment at KEK PQ . The analyzing power of the Collins FF was determined 
to be significantly non-zero, as can be seen in Figure l2*.121 This important measure- 
ment will provide vital input for factorized calculations, allowing pion production 
processes sensitive to the Collins mechanism to put constraints on the transversity 
distribution. 



Despite much progress within both the experimental and theoretical arenas in 
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Figure 2.12: Analyzing power of the Collins FF for pions as a function of z, measured 
by the BELLE experiment via two different methods. See pQ for details. 

recent years, no single, clear formalism has emerged in which to interpret the notable 
transverse spin effects that have been observed. A number of mechanisms remain 
as possible origins for the large transverse single-spin asymmetries, and a variety of 
further experimental measurements in different kinematic regions and sensitive to 
different partonic processes will be necessary in order to disentangle them. 
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Chapter 3 
QCD at RHIC 

RHIC was designed to study numerous aspects of QCD. The most flexible hadron 
collider in the world, it has so far produced gold-gold, copper-copper, deuteron- 
gold, and polarized proton-proton collisions at a variety of energies. Such a machine 
provides a rich environment for the study of QCD. The goal of colliding heavy ions 
at high energies is to create nuclear matter at extreme temperatures and densities, 
high enough that the quarks and gluons are (momentarily) not bound as hadrons but 
may co-exist rather as a quark-gluon plasma (QGP). Variations in collision species 
size and collision energy provide information on how the properties of the created 
matter are related to the initial conditions. Studies of small systems colliding with 
large nuclei (e.g. deuteron-gold) permit distinction between cold and hot nuclear 
effects. For a summary and review of what has been learned in the first few years of 
the RHIC heavy-ion program, see the evaluations from the four RHIC experiments 
published after the fourth year of running GHJ QUJ EI] . 

The nucleon structure program at RHIC, with unique access to high-energy 
polar ized-proton collisions, seeks to measure the helicity distributions of the par- 
tons within the proton, in particular gluon and sea-quark distributions, and to im- 
prove knowledge of the transverse spin structure of the proton. Through W boson 
production in eventual 500-GeV running, it will not only be possible and of inter- 
est to measure the flavor-separated helicity distributions of the sea quarks (Au(x), 
Ad(x)), but also the flavor-separated unpolarized pdf's (u(x), d(x)). As stated in 
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Section 12.1.41 a large difference in the u and d content of the (unpolarized) proton 
has been observed and is still not well understood. 

3.1 Cross section measurements and NLO pQCD 

The goal of this thesis is to demonstrate how proton structure can be investigated via 
proton-proton collisions at RHIC. It is essential to understand how well factorized 
pQCD can be used to describe and interpret the RHIC data. 

A number of polarization-averaged cross section measurements have been made 
at RHIC and compared to NLO pQCD calculations [TEJ E3 HE] • Comparisons 
of data to NLO pQCD calculations are shown in Figures 13.11 13.21 and 13.31 The 
most spectacular example of agreement between theory and data is seen in the 
mid-rapidity neutral pion measurement (Fig. 13. lj) |16j . which spans eight orders of 
magnitude and covers 1 < pr < 15 GeV/c. The data are compared to NLO pQCD 
calculations utilizing the CTEQ6M pdf 's |131j and two different FF sets, differing 
principally in the gluon-to-pion FF. While both calculations describe the data well 
down to what are perhaps surprisingly low values of pr (~ 1.5 GeV/c), the calcula- 
tion incorporating the FF set of Kniel, Kramer, and Potter (KKP) |113j is in better 
agreement with the data than that of Kretzer |114j . This is consistent with a larger 
gluon-to-pion FF. The bottom two panels in Fig. 13.11 indicate the sensitivity of the 
calculations to the choice of factorization and renormalization scales. The calcula- 
tions have been performed using equal factorization and renormalization scales of 
Pt, %Pt, and Pt/2. The KKP FF set in particular exhibits relatively little scale de- 
pendence. For comparison, the renormalization and factorization scale dependence 
of current DIS fixed-target experiments performing spin physics measurements is 
several times larger. An alternative theoretical calculation is compared to these 
neutral pion results in [7T] . 

There is similar agreement and similarly little scale dependence when studying 
the mid-rapidity production of inclusive charged hadrons, shown in Fig. 13.21 [15J. 
Here the calculations use the KKP FF set. As in the case of the neutral pions, the 
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Figure 3.1: Invariant cross section versus transverse momentum for mid-rapidity 
neutral pion production at PHENIX, taken from The data are compared to 
NLO pQCD calculations using two different gluon-to-pion fragmentation functions. 
See text for further details. 

bottom panel shows the difference between data and theory for three different scales 
of pt, 2pr, and pr/2. 

Even the forward production of neutral pions at RHIC, potentially susceptible 
to soft (non-perturbative) contributions, has been described by pQCD with relative 
success, as shown in Fig. 13.31 from the STAR collaboration j^j. The two calculations 
given in this figure utilize the KKP and Kretzer FF sets, and similar to the mid- 
rapidity data, the forward data, at least for p? ^1-7 GeV/c, are in better agreement 
with the KKP set. 
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Figure 3.2: Invariant cross section versus transverse momentum for mid-rapidity 
charged hadron production at PHENIX, compared to NLO pQCD, taken from |15j . 
See text for further details. 

The establishment of the ability of NLO pQCD to describe RHIC cross section 
data well and with little dependence on the choice of factorization and renormal- 
ization scales provides a solid foundation for using NLO pQCD to interpret in turn 
the polarized data at RHIC. 



3.2 Spin asymmetries in factorized QCD 



Generally, a spin asymmetry is the ratio of the difference to the sum of the spin- 
dependent cross sections for a particular process, given for example by Eq. 13. II in 
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Figure 3.3: Invariant cross section versus transverse momentum for forward neutral 
pion production at STAR, compared to NLO pQCD calculations using two different 
fragmentation functions, taken from See text for further details. 

the case of a double-spin asymmetry, with the arrow combinations representing 
same-spin and opposite-spin collisions which could be transverse or longitudinal. 



a 



TT 



a 



Tl 



(3.1) 



OrTt + o-TI 

The denominator is simply the total unpolarized cross section and as such is calcu- 
lable in factorized QCD as described above. The numerator is instead the difference 
of a convolution of spin-dependent pdf 's, a spin- dependent partonic hard scattering 
cross section, and spin-independent FF's. As in the polarization-averaged case, the 
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partonic hard scattering cross section is calculable directly from perturbative theory, 
while the pdf 's and FF's must be obtained from experimental measurements. The 
goal of the spin program at RHIC is to extract spin-dependent, or polarized, pdf's 
from asymmetry measurements. A particularly clean example is the extraction of 
the polarized gluon distribution, Ag, from the longitudinal double-spin asymmetry 
in direct photon production. At mid-rapidity and over a wide pt range at RHIC 
energies, about 75% of direct photon production comes from quark-gluon Compton 
scattering, q + g — > q + 7. Equation 13.21 gives the asymmetry in direct photon 
production from this process. 



Assuming the polarized quark distributions have already been well measured, e.g. in 
DIS experiments, it is relatively straightforward to extract Ag from Eq. 13.21 In 
practice, it is slightly more complicated, given that gluon Compton scattering is 
not the only process that contributes to direct photon production and the exact x 
values of the scattered partons are not known (see related discussion in Chapter EJ) • 
But generally, the longitudinal double-spin asymmetry in direct photon production 
at RHIC is expected to give a relatively clean measurement of Ag once enough 
statistics are available. 

The above example illustrates how factorized QCD can be used to extract po- 
larized pdf's from experimental spin asymmetries. In the future, once a variety of 
well-constrained measurements are available from RHIC, a global analysis of both 
DIS and RHIC results will be performed in order to obtain the spin-dependent pdf's 
that best describe all world data. 



A q L + r q+ \p T ) 



J2 q Ag(si) (g) Ag(x 2 ) Aq"+ 3 ^+7(g) 

E 9 qM ®g(x2) <8> (Ti+^i+^§) 



(3.2) 
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Chapter 4 
Experimental setup 

4.1 The Relativistic Heavy Ion Collider 

The Relativistic Heavy Ion Collider (RHIC) is located at Brookhaven National Lab- 
oratory on Long Island, New York. The RHIC storage ring is 3.83 km in circumfer- 
ence and is designed with six interaction points (IP's), at which beam collisions are 
possible. Up to 112 particle bunches per ring can be injected, in which case the time 
between bunch crossings at the IP's is 106 ns. The design luminosity for Au+Au 
collisions is 2 x 10 26 cm" 2 s~ 2 ; for p + p collisions it is 2 x 10 32 cm" 2 s -2 . The design 
polarization for proton beams is 70%. 

RHIC was built to collide heavy ions at a center-of-mass energy of up to 200 GeV 
per colliding nucleon pair and polarized protons at energies ranging from 50 to 
500 GeV. Collision of asymmetric species, i.e. different species in the two beams, 
is also possible due to independent rings with independent steering magnets. The 
first physics run took place in 2000, with Au+Au collisions at 130 GeV per nucleon. 
The following four running periods included Au+Au collisions at 200, 62.4, and 
19.6 GeV/nucleon, Cu+Cu collisions at 200, 62.4, and 22.4 GeV/nucleon, d+Au 
collisions at 200 GeV/nucleon, and polarized p + p collisions at 200 GeV. 

There were four major experiments developed for RHIC, three of which will 
continue to take data after 2005. There are two large experiments, STAR jHj and 
PHENIX [TJ], each of which have more than 500 collaborators, and two smaller 
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experiments, BRAHMS |5J and PHOBOS jSU], having fewer than 100 collaborators 
each. PHENIX, STAR, and BRAHMS all have a spin-physics program as well 
as a heavy-ion program, while PHOBOS only studies heavy-ion physics. PHOBOS 
finished taking data in 2005. The four experiments were designed with some overlap 
and some complementarity in the physics processes they could measure. In this way 
it is frequently possible for one experiment to cross-check the results of another, 
yet each experiment has its own area of specialization. The PHENIX experiment, 
through which the measurement for this thesis was made, is described in detail in 
Section 14.31 below. 

4.2 RHIC as a polarized p+p collider 

RHIC is the first and only polarized proton-proton collider in the world. Figure l4~T1 
shows a diagram of the RHIC complex including all equipment relevant for polarized 
proton beams. A number of technological developments and advances have made it 
possible to create a high-current polarized source, maintain the beam polarization 
throughout acceleration and storage, and obtain accurate measurements of the de- 
gree of beam polarization at various stages from the source to full-energy beams in 
RHIC. For an overview of RHIC as a polarized-proton collider, refer to [2*K] . 

4.2.1 RHIC- AGS complex 

In the case of polarized-proton running at RHIC, a pulse of polarized H~ ions from 
the source (see Section l4.2.2j) is accelerated to 200 MeV in the Linac, then stripped 
of its electrons as it is injected and captured as a single bunch of polarized protons 
in the Booster, which accelerates the protons to 1.5 GeV. The bunch of polarized 
protons is then transferred to the Alternating Gradient Synchrotron (AGS) and 
accelerated to 24 GeV before injection into RHIC. Each bunch is accelerated in the 
AGS and injected into RHIC independently, with the two RHIC rings being filled 
one bunch at a time. The direction of the spin vector is selected for each bunch 
separately. 
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Figure 4.1: The RHIC- AGS complex as a polarized proton facility. 
4.2.2 Polarized source 

Polarized proton injection uses an optically-pumped polarized H~ ion source (OP- 
PIS) constructed at TRIUMF in Canada from an OPPIS source previously used at 
KEK in Japan |148j . Polarization of 85% has been reached. In order to achieve a 
polarized proton beam, there are several steps. Unpolarized protons first pick up 
polarized electrons from an optically pumped rubidium vapor in a high magnetic 
field, forming hydrogen atoms with polarized electrons. The electron polarization 
is then transferred to the nucleus via static magnetic fields, creating nuclear-spin- 
polarized atoms. These atoms pick up a second unpolarized electron in a sodium 
vapor negative-ionizer cell. A pulse of polarized H~ ions is then accelerated to 
200 MeV in the Linac and injected as a single bunch of polarized protons in the 
Booster. 
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4.2.3 Siberian snakes 

The precession of the spin vector of a proton in a planar circular accelerator ring is 
given by the Thomas-BMT equation [T4HI57] (Eq.lO|). 



in which e is the proton charge, M is the proton mass, 7 = E/M is the relativistic 
boost, G = 1.7928 is the anomalous magnetic moment of the proton, B± indicates 
the magnetic field perpendicular to the plane of proton motion, typically the vertical 
guide field, and B» is the longitudinal field. P here is in the frame of the proton. 
Note that for highly energetic protons (7 large), the B± term dominates. 

Comparing the spin evolution equation to the Lorentz force equation of motion 
for a particle orbit in a magnetic field (Eq. 14.2)1 . 



it becomes evident that for highly energetic protons or in a purely vertical field, the 
spin precesses a factor of G7 faster than the orbital motion, meaning that the spin 
precesses G7 times in a single revolution around the RHIC ring. This number is 
referred to as the spin tune, u sp . 

There are two principal types of depolarizing resonances that may be encountered 
during acceleration, imperfection and intrinsic resonances. Imperfection resonances 
refer to those due to errors in magnet currents or alignments; intrinsic resonances are 
due to the (desired) beam focusing fields. A depolarizing resonance is encountered 
whenever the spin precession frequency is equal to the frequency with which a de- 
polarizing field is crossed, leading to additive depolarization effects. Imperfect fields 
exist with a more or less random distribution around each ring. As such, resonances 
occur when v sp = G7 = n, an integer, and the spin vector is at the same phase in 
its precession every time it encounters the same imperfect field. Intrinsic resonances 
occur when u sp = G7 = kP ± i/ y , where k is an integer, P is the superperiodicity, or 
degree of regularity in the focusing-defocusing lattice, and u y is the vertical betatron 
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tune, the number of oscillations around the stable beam orbit per beam revolution, 
in the vertical plane. (The z-axis is taken to be in direction of proton motion, in 
the frame of the proton.) 

Siberian snakes a series of spin- rotating dipoles, so named because of the 
beam trajectory through the magnets and the fact that they were developed at 
Novosibirsk, are used to overcome both imperfection and intrinsic resonances in 
RHIC. There are two snakes installed in each RHIC ring at diametrically opposite 
points along the rings. The two snakes in each ring rotate the spin vector 180° about 
perpendicular horizontal axes, without perturbing the stable spin direction and with 
only local distortion of the beam orbit. In this way, all additive depolarizing effects 
from resonances are avoided. 

In the AGS there is not enough space to permit a full snake. Only a partial snake 
[134 , which rotates the spin vector by less than 180°, is possible, and complete re- 
versal of the spin direction only occurs over the course of multiple revolutions. A 
partial snake of less than 10% is sufficient to overcome imperfection resonances but 
not intrinsic resonances in the AGS; a partial snake of 20% or more is expected to 
overcome both types of resonances. In the 2001-02 run, a solenoidal "warm" (non- 
superconducting) snake was available in the AGS. This magnet was a 5% snake; 
therefore, it rotated the spin direction by 9° each revolution. In order to handle 
intrinsic resonances, the technique used was to artificially enhance the resonances 
such that they were tuned to produce a complete spin flip each time one is encoun- 
tered, rather than depolarization. An AC dipole magnet was used to achieve this. 
In 2004 the solenoidal snake was upgraded to a warm helical snake, which aided in 
decoupling the snake's effects on the x and y motion of the beam. 

In 2005 a superconducting "cold" snake was installed in the AGS and partially 
commissioned; this final piece of equipment should allow achievement of the RHIC 
design polarization of 70%. Superconducting magnets allow a larger magnetic in- 
duction, thus a larger j B ■ dl and greater spin rotation per revolution. The cold 
snake is a 25% snake, rotating the spin vector 45° each time. In 2006 the plan 
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is to use both the warm and cold snakes in conjunction, rotating the spin by 5% 
and 15% respectively, achieving the minimum rotation of 20% required to handle 
all resonances in the AGS. A description of the technique of using multiple partial 
snakes in conjunction can be found in |135| . 

The first polarized-proton run in 2001-2002 achieved an average beam polariza- 
tion in RHIC of ~ 15%. Improvements in both available hardware and machine 
understanding led to measured polarization values in excess of 50% in 2005. In ad- 
dition, significant polarization (~ 30%) at a center-of-mass energy of 410 GeV was 
achieved in 2005. Further high-energy commissioning in preparation for eventual 
running at y/s = 500 GeV will be done in 2006. 

4.2.4 Polarimeters 

For RHIC to provide full-energy polarized beams, the polarization must be mea- 
surable at various stages of acceleration in order to identify and address possible 
origins of depolarization at each step. Only RHIC polarimetry will be discussed 
here and not the various techniques used to measure the beam polarization further 
upstream in the process of acceleration. The RHIC polarimetry fulfills a three-fold 
purpose: (1) beam polarization measurements to provide feedback for the acceler- 
ator physicists; (2) beam polarization measurements to provide polarization values 
as input for the various experiments; and (3) experimental study of polarized elastic 
scattering. 

There are two types of polarimeters installed in RHIC, designed to measure the 
beam polarization in the vertical direction. The proton-carbon (pC) polarimeter 
takes advantage of a known analyzing power, A p ^ m 0.01, in the elastic scattering 
of polarized protons with carbon atoms (p^ + C — > p^ + C) in order to measure 
the beam polarization. This analyzing power originates from interference between 
electromagnetic and hadronic elastic scattering amplitudes, which are finite due to 
the anomalous magnetic moment of the proton. Thus the term "Coulomb-nuclear 
interference" (CNI) is used to describe the process. Refer for example to (ZOj for a 
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discussion of CNI in hadronic reactions at high energies. For the pC polarimeter at 
RHIC, the left-right (azimuthal) asymmetry of the recoil carbon atoms is measured 
using an array of silicon detectors, allowing calculation of the beam polarization as 
in Eq. |OJ 



where N L (N^) indicates the number of recoil carbons scattering to the left (right) 
of the proton beam from bunches polarized up (down). (Compare to Eq. 15.81 in 
Section EXU) 

The analyzing power in the process now utilized by the pC polarimeter to mea- 
sure the RHIC beam polarization was initially measured by AGS experiment E950 
|26j . However, the AGS measurement was only made to ±30%, at a beam en- 
ergy of 22 GeV, and further calibration is necessary to understand the RHIC beam 
polarization at 100 GeV. 

Calibration of the pC polarimeter to within an absolute beam polarization of 
~ 5% can be provided by a polarized hydrogen-jet-target polarimeter |149j . Such a 
device was commissioned at RHIC in 2004. The measurement is in some ways simi- 
lar to that of the pC polarimeter. The left-right transverse single-spin asymmetry of 
elastically scattered protons due to the CNI process is measured using silicon strip 
detectors. However, the analyzing power for this process of elastic proton-proton 
scattering is not known. Instead, the hydrogen-jet-target polarization, > 90%, is 
known to better than 2% in absolute polarization. The beam polarization can be 
obtained by measuring first the asymmetry due to the polarized target (stgt), aver- 
aging over the spin states of the beam, and then similarly measuring the asymmetry 
due to the polarized beam (e^am); averaging over spin states of the target, which are 
varied in time. The same (unknown) analyzing power, An, applies to both cases, 



beam 




(4.3) 
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and the beam polarization, Pbeam, can be calculated as in Eq. 14.41 



A N 




beam 



beam tgt 



£ tgt 



(4.4) 



Once the beam polarization is determined, Eq. 14 .HI can be rearranged and used in 
turn to measure the analyzing power of the process. Results on the measurement of 
the analyzing power of the jet polarimeter can be found in |127j . 

Low rates for this process mean that measurements performed over a long time 
are necessary. Because of this, the pC polarimeter, which can make measurements 
in less than ten seconds and provide immediate information on the stability or decay 
of the beam polarization from a few data points taken over the several hours of a 
fill, remains essential. Results obtained from the two different polarimeters can be 
compared offline, and in this way the pC polarimeter can be calibrated. 

In order to check the absolute sign of the spin direction for each bunch, it is pos- 
sible to measure bunch-by-bunch asymmetries, making left-right asymmetry mea- 
surements from events coming from only a single bunch number, and then take 
advantage of the fact that the sign of the CNI asymmetry is known from theory. 

Note that if bunch-by-bunch polarimeter information were not available, as was 
the case in 2002, and a few individual bunches had the incorrect spin direction 
assigned to them, this would make the beam polarization appear to be less than 
its true value. The experiments in turn would utilize these same, incorrect spin 
assignments, diluting any raw asymmetry in particle production they may observe. 
However, the larger beam polarization correction would exactly compensate for the 
diluted raw asymmetry, and a true physics asymmetry would still result. If instead 
all bunches were assigned the incorrect spin direction, the measured CNI asymmetry 
would have the opposite sign with respect to theoretical predictions. 

4.2.5 Spin rotators 

The naturally stable spin direction through acceleration and storage is transverse 
to the proton's momentum, in the vertical direction. Spin rotator dipole magnets, 
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commissioned in 2003, have been used to achieve both radial and longitudinal spin 
|121j . The spin vector is rotated away from vertical immediately before the collision 
point and then back to vertical immediately afterwards. The rotators are located 
outside the interaction regions of the PHENIX and STAR experiments, giving both 
experiments the ability to choose independently whether they want longitudinally 
or transversely polarized collisions. The BRAHMS experiment, having no spin ro- 
tators available, focuses on transverse spin measurements. The local nature of the 
spin rotator magnets means that the STAR and PHENIX experiments must each 
have their own way of checking the direction of the spin vector at their respective 
interaction regions. In PHENIX this ability is provided by the local polarimeter, 
described in Section 14.3.31 

4.2.6 Spin flippers 

In addition, RHIC is equipped with a "spin flipper," an AC dipole with a radial 
magnetic field installed in an area common to both beams, which flips the spin 
180° bunch by bunch Spin flipping involves detuning one of the two snakes to 
alter the spin tune and then sweeping the spin flipper frequency of the AC dipole 
through the resonance to effect a controlled spin reversal. The ability to flip all the 
spin vectors bunch by bunch is important in order to reduce potential systematic 
errors from any correlations that may exist between a bunch and its spin direction. 
An example of such possible correlations would be a systematically larger bunch 
length for bunches of a particular spin sign, which could lead to subtle differences in 
the effective luminosity seen by the luminosity detectors and the final-state-particle 
detectors (e.g. the PHENIX central arms, see Section l4.3.2j) . Frequent or carefully 
timed use of the spin flipper throughout a machine fill could also facilitate balancing 
the relative luminosity at all interaction regions; any bunch crossings later discarded 
for any reason from offline analysis would represent the various spin configurations 
in an approximately balanced way, and the relative luminosity would not be greatly 
affected. 
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Preliminary commissioning of the spin flippers was performed in 2002 [S3] but 
was not continued in subsequent years. Final commissioning and the beginning of 
regular usage is planned for 2006 or 2007. 

4.3 The PHENIX experiment and detector 

The PHENIX collaboration, one of the two large collaborations at RHIC, is com- 
prised of more than 500 people from over 60 institutions around the world. The 
PHENIX experiment was designed to specialize in the measurement of photons, 
electrons, and muons as well as high-py probes in general over a limited accep- 
tance, with good mass resolution and particle identification capabilities. It has an 
extremely high rate capability as well as sophisticated trigger systems, allowing mea- 
surement of rare processes. PHENIX produced a wealth of physics results from the 
first years of RHIC running and anticipates the ability to continue taking advantage 
of improvements in machine luminosity to further develop its program to measure 
rare probes. 

The PHENIX detector consists of two mid-rapidity (|^| < 0.35, i] = — In tan |) 
spectrometer arms, primarily for identifying and tracking charged particles as well 
as measuring electromagnetic probes, forward spectrometer arms for identifying and 
tracking muons, and interaction detectors. An overview of the PHENIX detector is 
given in [T2"] . 

4.3.1 Interaction and vertex detectors 

PHENIX has two detectors used to determine when a collision occurs and to mea- 
sure the position of the event vertex in the beam direction. The beam-beam coun- 
ters (BBC's) are quartz Cerenkov detectors, located at ±1.44 m from the nominal 
interaction point, which measure charged particles over 2-7T and cover a pseudora- 
pidity range of 3.0 < \q\ < 3.9. The zero-degree calorimeters (ZDC's) are hadronic 
calorimeters sensitive primarily to forward neutrons. They also cover 2tt in az- 
imuth, are located at ±18 m, and extend over 4.7 < \rj\ < 5.6. The ZDC's were 
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developed for RHIC rather than specifically for PHENIX and are used by all four 
RHIC experiments. 

A minimum-bias (MB) trigger occurs when there is a minimum of one photo- 
multiplier tube (PMT) fired in each of the two BBC's. In 200-GeV proton-proton 
collisions, the BBC's see approximately 50% of the total inelastic p + p cross section 
and provide a vertex resolution of ~ 2 cm. The ZDC's are sensitive to a much 
smaller fraction of the interactions and are not used as part of the MB trigger in 
p + p. The primary roles of the ZDC's in proton-proton collisions are in local po- 
larimetry (see Section 14.3.3)) and for systematic studies of the relative luminosity 
between bunch crossings with different spin configurations (see Appendix E}. In 
gold-gold collisions both the BBC's and ZDC's observe a significant fraction of the 
inelastic cross section (the BBC's observe 92%), and correlations in the measured 
ZDC energy and BBC charge are used to determine the centrality (degree of overlap) 
of the heavy-ion collisions. 

In the 2001-02 p + p run at RHIC, when the data for the present work were 
taken, there was an additional interaction detector installed, the normalization trig- 
ger counter (NTC). The NTC consisted of two scintillation counters located at 
±40 cm from the center of the interaction diamond and covering the full azimuth. 
This detector was installed to improve cross section measurements in p+p collisions; 
it was sensitive to ~ 85% of the total inelastic cross section. For further discussion 
of the PHENIX BBC's and NTC, see j2H]; for a more detailed description of the 
RHIC ZDC's, see [Hj. 

4.3.2 PHENIX central arms 

The PHENIX central arms cover a pseudorapidity range of \rj\ < 0.35 and two 90- 
degree intervals in azimuth, offset 33.75 degrees from vertical. See Figure l4~2l for 
a schematic diagram of the PHENIX central arms as instrumented in the 2001-02 
RHIC run. 
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PHENIX Detector - Second Year Physics Run 




West Beam View East 

Figure 4.2: PHENIX central arms in 2001-02. 
4.3.2.1 Tracking detectors 

The PHENIX central tracking system is comprised of a drift chamber (DC) and pad 
chambers (PC's). The DC is located at r = 2.0 — 2.4 m from the beam pipe and 
covers 2 m along the beam direction. It measures the trajectories of charged particles 
in the r — <ft direction, providing the px of the particle as well as <ft information, 
necessary for invariant-mass reconstruction of particle pairs. The innermost plane 
of the PC is located immediately outside the DC; the outermost plane is immediately 
in front of the electromagnetic calorimeter (EMCal, see below). There is a middle 
PC plane located behind the ring-imaging Cerenkov detector (see below) in the West 
arm only. The PC's provide three-dimensional spatial point information for pattern 
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recognition. The outermost plane of the PC is particularly effective in providing a 
charged veto for particles that shower in the EMCal. 

The tracking detectors are placed outside the magnetic field; there is only a 
residual field of < 0.6 kG present in the DC. With no inner tracking available, 
tracks are assumed to point back to the event vertex, provided by the BBC's. Thus 
charged particles that do not originate at the vertex have incorrectly reconstructed 
momentum, leading to low-momentum, long-lived particle decays and conversion 
electrons as sources of background in high-p-^ charged particle analysis. The tracking 
detectors are described in more detail in [TT] . 

4.3.2.2 Ring-imaging Cerenkov detector 

A ring-imaging Cerenkov detector (RICH) is placed between the layers of tracking 
subsystems. Its main purpose is electron-pion discrimination. It has been filled 
with CO2 gas in the first five years of RHIC running. The momentum threshold 
for production of Cerenkov radiation by charged pions in CO2 is 4.7 GeV/c, while 
for electrons it is 0.017 GeV/c. In the pr range between these two values, there is 
excellent e/ir separation. Further discussion of the RICH can be found in [2"3*] . 

4.3.2.3 Electromagnetic calorimeter 

The EMCal is designed to measure the energy and position of electrons and photons. 
It also contributes significantly to particle identification through energy-position 
measurements as well as timing information. It is the outermost subsystem in the 
central arms, located at ~ 5 m in r. The EMCal in the West arm is divided 
into four sectors of a lead-scintillator (PbSc) sampling calorimeter; in the East arm 
there are two sectors of PbSc and two of a lead-glass (PbGl) Cerenkov calorimeter. 
The two kinds of detectors rely on significantly different physics processes and as 
such can be useful in making systematic comparisons. The PbGl excels in energy 
measurements and is less sensitive to hadrons, while the PbSc has better timing 
resolution. The nominal energy resolution of the PbSc is 8.1% /WE (GeV) © 2.1%; 
the PbGl has a nominal energy resolution of 6%/ a/ E (GeV). The intrinsic timing 
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resolution is ~ 200 ps for the PbSc and ~ 300 ps for the PbGl. While designed as 
an electromagnetic calorimeter, the EMCal is approximately 1 nuclear interaction 
length and offers some sensitivity to hadrons. Detailed description of the EMCal is 
available in |4*T] . 

4.3.3 Local polarimeter 

In order to provide local polarimetry, a shower maximum detector (SMD) was added 
to the PHENIX ZDC's. The ZDC/SMD combination measures both the energy and 
position of forward neutrons. A large (~ —11%) transverse single-spin azimuthal 
asymmetry for neutrons with < 0.5 GeV/c has been observed [SU]- This asym- 
metry can be exploited to measure the degree to which the beam polarization is 
vertically transverse, radially transverse, or longitudinal. Maxima and minima of 
particle production are expected at right angles in to the spin vector. For longi- 
tudinal polarization, no azimuthal asymmetry is expected; a non-zero longitudinal 
single-spin asymmetry would violate parity and is forbidden in particle production 
via strong and electromagnetic processes. Figure 14.31 shows the observed asymme- 
try as a function of <fi for the cases of vertical and longitudinal beam polarization. 
The expected maxima and minima at ±| can be seen for vertical polarization, and 
the asymmetry is approximately consistent with zero in the case of longitudinal 
polarization. 

4.3.4 Muon arms 

The forward spectrometer arms in PHENIX are designed to identify and measure 
prompt and decay muons. They cover 2tt in azimuth and a pseudorapidity range of 
1.2 < < 2.4. A high-resolution tracker in a radial magnetic field is followed by 
alternating layers of absorber and low-resolution tracking for muon identification. 
The tracker magnet also serves as a hadron absorber. The design is to reject pions 
and kaons at a level of 10~ 3 . In Figure fOl the muon arms are shown as instrumented 
for the 2001-02 RHIC run. Note that only the south arm was available for data 
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Figure 4.3: PHENIX local polarimeter measurement. The raw asymmetry in neu- 
tron production divided by the degree of beam polarization is shown vs. azimuthal 
angle. The solid points are with the spin rotators off (vertical polarization); the 
open points are with the spin rotators on (longitudinal polarization). 



taking in this year; the north arm was completed in 2003. See 
information regarding the muon arms. 



for further 



4.3.5 Detectors and triggers used in the analysis for this 
thesis 

To measure the energy and position of photons from neutral pion decays, the present 
analysis utilized both the PbSc and PbGl electromagnetic calorimeters. The cen- 
tral arm tracking detectors were used to eliminate EMCal clusters associated with 
charged tracks. The BBC's were used to provide a MB trigger. The EMCal-RICH 
trigger (ERT, described below) was used to select events with high-p T photons. 

4.3.5.1 EMCal-RICH trigger 

The hardware- level rare-event triggers implemented in PHENIX include an EMCal- 
RICH trigger (ERT). The ERT can be used as a high-energy photon trigger or in con- 
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PHENIX Detector - Second Year Physics Run 




□ Installed □ Active 

Figure 4.4: PHENIX muon arms in 2001-02. 

junction with the RICH as an electron trigger. The EMCal towers are grouped into 
trigger tiles, over which the deposited energy is summed, and an event is recorded 
if the total energy in a tile is above a preset threshold. The tiles are configured to 
be non-overlapping sets of 2 x 2-tower regions and overlapping sets of 4 x 4-tower 
regions. Up to three different energy thresholds may be set for the 4 x 4-tower tile 
size. The electron trigger utilizes the 2 x 2-tower tile size. 

4.3.5.2 EMCal calibration 

There was an initial series of EMCal calibrations performed, both before detector 
installation and in situ in the PHENIX interaction region, descriptions of which can 
be found in jH] . During regular data taking, a pulsed UV laser is utilized to provide 
EMCal calibration information. In offline analysis, three major techniques used to 
calibrate the EMCal energy are examining the location of the 7r° mass peak position, 
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the energy-momentum ratio for electrons, and the location of the minimum-ionizing 
peak. 
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Chapter 5 

Data and analysis 

5.1 Overview 

The left-right transverse single-spin asymmetry, A at, of mid-rapidity neutral pions 
was measured from the data taken during the first polarized proton run at RHIC, 
in late 2001 and early 2002. To make a single-spin measurement with two polarized 
beams, the spin states of only one beam at a time were taken into account, averaging 
over the spin states of the other. Neutral pions were reconstructed via their decay 
to two photons. The general procedure was to obtain the spin-dependent ir° yields 
for each machine fill, calculate the raw (uncorrected) asymmetries for each fill as 
a function of pt, make fill-by-fill polarization corrections, and then average the 
asymmetries over all fills. The results were corrected for estimated contributions to 
the asymmetry from the background under the 7r° invariant-mass peak by measuring 
and subtracting the asymmetry of the background immediately around the peak in 
mass. Various studies were subsequently performed as checks. The final results of 
this analysis have been published in Physical Review Letters ^H] • 

5.2 Data selection and quality 

Prior to this spin-asymmetry analysis, a polarization-averaged cross section mea- 
surement was made for neutral pion production at PHENIX and published in jl6j . 
The analysis completed for this thesis started with the data sample utilized in the 
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cross section analysis, and then additional quality cuts relevant to a spin-dependent 
analysis were made. 

The data used were from ERT 2 x 2-triggered events. The 2x2 energy threshold 
was ~ 0.8 GeV in the 2001-02 run. The triggered sample had much better statis- 
tics at higher transverse momenta than the MB sample, but the MB sample was 
also analyzed for comparison purposes. 18.7 million triggered events were analyzed, 
corresponding to approximately 880 million sampled MB events. The detector sub- 
systems involved were the BBC and the EMCal, including both the lead scintillator 
(PbSc) and lead glass (PbGl). Additionally, the drift chamber and pad chambers 
were used to veto clusters in the EMCal produced by charged particles. Basic checks 
on the quality of output from all involved subsystems were performed. 

5.2.1 Fill and run selection 

In making 7r° —+7 + 7 spin-asymmetry or cross section measurements, it is essential 
to ascertain that the EMCal was working properly when the analyzed runs were 
taken. The uppermost sector of PbSc in the West arm had been omitted from the 
cross section analysis due to problems with electronics noise during data taking. A 
number of runs had been eliminated from that analysis because of EMCal towers 
with unusually high or low numbers of hits that were not among the known and 
understood hot or dead EMCal towers. In addition, a sequence of runs had been 
omitted because of the ERT. The energy threshold for the ERT 2x2 trigger had 
been adjusted in the early part of the data-taking period, and only runs taken after 
the threshold was stable at ~ 0.8 GeV were included. Runs from a total of twenty 
machine fills were included in the final cross section analysis. 

For the present analysis, six of the twenty fills used in the cross section analysis 
were removed. Two were removed because no polarization measurement was avail- 
able. Three additional fills were removed because a number of bad runs from these 
fills were found in a study investigating the stability of the MB trigger. One fill was 
removed because it had an unusual spin pattern. 
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5.2.2 Event and crossing selection 

An offline BBC event vertex within ±30 cm of the nominal interaction point was 
required for all events. The acceptance of the central arms is approximately constant 
for collisions taking place in this region. The BBC vertex resolution was ~ 2 cm. 
The online vertex cut for the MB trigger in the 2001-02 run was ±75 cm. 

Spin-sorting the events requires keeping track of the direction of the spin vector 
for the polarized bunch in the bunch crossing that produced the event (simply the 
term crossing will be used henceforth). Note that the same pairs of bunches in the 
two rings collide at the same interaction points each time, but the same pairs of 
bunches do not necessarily collide at the different interaction points, i.e. at different 
experiments. Data quality was checked on a crossing-by-crossing basis. For the 
entire 2001-02 p ± p data set, four crossings out of a nominal 60 were consistently 
removed from this analysis. These four crossings had unusually low luminosities 
because they were regularly affected by either injection or steering activities in the 
ring. In addition, the ten crossings in which only one beam had filled bunches 
while the other had a five-crossing beam-abort gap were removed from all fills. Any 
events occurring in these crossings were beam-background collisions rather than 
beam-beam collisions. 

There were two bunch- by-bunch spin patterns utilized in the 2001-02 run, one 
in the initial period of data taking and the other in the latter period. The second 
pattern included bunches with zero polarization so that they would be available for 
systematic checks, in particular for the pC polarimeter. For fills with the second po- 
larization pattern, one additional crossing was eliminated from this analysis because 
it had zero polarization. The other zero-polarization crossings corresponded to two 
of the low-luminosity bunches that had already been eliminated. A schematic illus- 
tration of the first spin pattern is given in Figure 15.11 Note that the spin patterns 
are selected to provide approximately equal numbers of same-spin and opposite-spin 
crossings for each experiment. They are also chosen so that the spin combinations 
are different every crossing. This rapid change in spin combinations, occurring on 
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Figure 5.1: Schematic illustration of the first spin pattern, with no zero-polarization 
crossings. In one beam, alternate bunches have opposite spin directions; in the other 
beam, the spin direction changes every two bunches. 

the order of every hundred nanoseconds, greatly reduces potential time-dependent 
systematic uncertainties. 

Crossing-by-crossing luminosity measurements from the pC polarimeter provided 
a fill-dependent bad crossing list. A total of four individual crossings from two fills 
were discarded from the analysis based on this information. 



5.3 EMCal-RICH trigger 

All data used for the final asymmetry results were from events accepted by the 
2 x 2-tower tile trigger for high-p^ (high energy at mid-rapidity) photons. The 2x2 
trigger had an average rejection factor of 47, i.e. only accepted on average one in 47 
MB events. It had a 78% efficiency for neutral pions above px ~ 3.5 GeV/c, as can 
be seen in Fig. 15.21 

The effects of the trigger on the mean px of the photon pairs falling under the ir° 
peak can be seen in Table IB~T1 The trigger raised the mean px in the 1-2 GeV/c bin 
significantly because it was still well below its maximum efficiency in this transverse 
momentum range. The trigger had little effect on the mean px of higher px bins. 
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Figure 5.2: a) The 2x2 trigger efficiency for neutral pions, as a function of pion 
transverse momentum. The dashed line shows a Monte Carlo simulation based on 
trigger tile efficiencies, and the solid line indicates an upper limit on the 7r° efficiency 
based on the number of active trigger tiles, b) The fraction of the 7r° yield satisfying 
the MB trigger condition. The solid line is the fit of the data to a constant. The 
figure is taken from |16j . 
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1-2 


1.27 


1.38 


2-3 


2.32 


2.33 


3-4 


3.33 


3.35 


4-5 


4.39 


4.36 



Table 5.1: Mean pt of photon pairs under the ir mass peak for MB and 2x2- 
triggered data. 
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V t (GeV/c) 


7r° yield 


Peak width (MeV/c 2 ) 


1-2 


658k 


13.2 


2-3 


143k 


11.2 


3-4 


22k 


10.4 


4-5 


4k 


10.6 



Table 5.2: 7r yields obtained after background subtraction; la tt peak widths from 
a Gaussian fit. 

5.4 Reduction of background 

In order to understand the background contribution to the 7r° mass peak and obtain 
7T° yields, the invariant-mass spectrum for photon pairs was fitted. The mass peak 
was fitted to a Gaussian and the combinatorial background to a second-degree poly- 
nomial. The 7T° yields per p? bin, given in Table 15.21 were obtained by subtracting 
the background from the total number of pairs in the peak. An example fitted and 
subtracted invariant-mass spectrum for 1 < px < 2 GeV/c is shown in Figure I5~31 
Mass bins of 10 MeV/c 2 were used in this analysis. 

The la peak widths from a Gaussian fit to the 120-160 MeV/c 2 mass region are 
shown in Table 15.21 The transverse momentum spectrum for the pairs falling under 
the 7T° mass peak can be seen in Figure 15.41 

Several cuts were made in order to reduce the background in the photon pair 
sample. Mismatched true photons, coming from two different particles, were the 
main source of background. Other sources of background included electrons, hadrons 
that deposited energy in the EMCal, or secondary particles not coming from the 
event vertex, all of which could lead to false combinatorial pairs whose mass fell 
under the 7r° mass peak, either in combination with each other or with true photons 
from neutral pions. The cuts were: 

• Minimum energy cut of 0.1 GeV in the PbSc, 0.2 GeV in the PbGl. This cut 
effectively eliminated pairs in which one photon carried nearly all the energy 
and the other very little. The same minimum energy cut was used for the ir° 
An analysis, the results of which were published in [TTj . 
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Figure 5.3: Top panel: Invariant-mass spectrum for 1 < pt < 2 GeV/c photon pairs. 
Middle panel: Fitted spectrum. Bottom panel: Subtracted spectrum. 

• Charged veto cut. All EMCal clusters within a 10-cm radius of a projected 
charged-track position onto the EMCal were excluded. 

• Shower shape cut in the PbSc to select clusters displaying the expected shape 
for energy deposits from photon hits. The analogous information for the PbGl 
was not available. 



While a time-of-flight (TOF) cut could have potentially offered additional hadron- 
photon discrimination, the EMCal TOF was not well calibrated in the 2001-02 data 
set, so no timing cut was performed. 
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p T (GeV/c) 


Background fraction (%) 




Before cuts 


After cuts 


1-2 


58 


34 


2-3 


23 


12 


3-4 


12 


6 


4-5 


9 


5 



Table 5.3: Reduction of background contribution to the tt mass peak, taken as 
120-160 MeV/c 2 , before and after background-removal cuts. 



The background fraction, r, was obtained by taking the ratio of the fitted back- 
ground to the total number of pairs falling within the mass peak. This fraction 
before and after the extra cuts were performed is shown in Table 15.31 Significant 
reduction was achieved in the lowest two pt bins. 
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5.5 Asymmetry calculation 
5.5.1 Overview 

The two counter-circulating RHIC beams are frequently referred to as "blue" and 
"yellow," named after the colored stripes painted on the respective magnet systems. 
The blue beam orbits clockwise, the yellow beam counter-clockwise. Both beams 
are typically polarized, as they were in the 2001-02 data-taking period. In order 
to make a single-spin asymmetry measurement, the spin direction of the bunches 
in only one beam was considered at a time, averaging over the spin direction of 
the bunches in the other. Results from the blue and yellow beams were obtained 
separately and subsequently combined. 

The formula used to calculate asymmetry values is given in Eq. 15.11 

i i - 

N ~ p hc ^(\cos^\)m + nm { ■ j 

in which Pbeam is the beam polarization, ^ c ^ sip ^ is an azimuthal acceptance correc- 
tion factor (see below), (N^) is the neutral pion yield from crossings with the 
polarized bunch spin up (down), and 1Z = £} j £} is the relative luminosity between 
crossings having the polarized bunch with spin up versus down. As An is a left- 
right asymmetry, Eq. 15.11 must be used separately for the two detector arms. As 
given, it applies to yields to the left of the polarized beam; an overall minus sign is 
required for yields observed to the right of the polarized beam. Asymmetry results 
for the left and right detector arms were obtained separately and then combined. 
The correction factors and relative luminosity are discussed further below. 

The beam polarization varies fill by fill. Thus the asymmetry is determined for 
every fill, then averaged over all fills. An example of fill-by-fill asymmetries is given 
in Figure 15 .5| this figure also gives an indication of the fill-by-fill stability of the 
measured asymmetry. Large observed variation among fills could indicate system- 
atic errors. (Note that all uncertainties given in figures and tables are statistical 
unless stated otherwise.) Table E31 shows results obtained from Eq. 15.11 for the two 
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Figure 5.5: Fill-by-fill asymmetry results for the two beams. 

detector arms for triggered events. Note that the uncertainties are slightly smaller 
for the yellow beam due to higher average polarization. Combined results for both 
detector arms are shown in Figure I5.6[ and the values for both arms and beams 
combined are given in Table I5.5I All results given in this section are before correc- 
tion to the asymmetries for the asymmetry of the background, which is described 
in Section IB. 5. 51 



5.5.2 Determination of relative luminosity 

The relative luminosity between crossings with spin-up bunches and crossings with 
spin-down bunches for the polarized beam in consideration was obtained from the 
number of MB events recorded by the BBC. A typical value in this analysis was TZ = 
£ T /C l = 1.09 for the yellow beam, K = 0.92 for the blue beam. Error in the relative 
luminosity as measured in this case for azimuthal transverse single-spin asymmetries 
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vt (GeV/c) 


Beam 


Left 


Right 






Am 


(T A 


Am 


(T A 


1-2 


Blue 


-0.008 


0.015 


-0.012 


0.019 




Yellow 


-0.005 


0.015 


0.003 


0.011 


2-3 


Blue 


0.028 


0.035 


-0.066 


0.039 




Yellow 


-0.016 


0.030 


-0.003 


0.028 


3-4 


Blue 


-0.101 


0.094 


0.106 


0.099 




Yellow 


0.033 


0.077 


-0.092 


0.073 


4-5 


Blue 


-0.02 


0.23 


0.16 


0.22 




Yellow 


0.01 


0.17 


0.14 


0.18 



Table 5.4: Asymmetry results for particles observed in the left detector arm (west 
arm for blue, east arm for yellow) and right detector arm (east for blue, west for 
yellow) . 



Vt (GeV/c) 


A N 


a A N 


1-2 


-0.006 


0.008 


2-3 


-0.014 


0.017 


3-4 


-0.013 


0.043 


4-5 


0.070 


0.101 



Table 5.5: Combined asymmetry results for the two beams and particle production 
on the two sides of the polarized beam. 

could potentially come from an azimuthal dependence of the BBC efficiencies. The 
error would be proportional to the physics asymmetry of the particles hitting the 
BBC (measured to be < 1%) times the difference in the efficiency of the left and 
right halves of the detector. There is no current measurement of this value, but it is 
expected to be small, also at the level of a few times 10~ 2 or less. This would lead 
to a potential systematic error on 1Z on the order of a few times 10~ 4 , but this is 
only a rough estimate. Rather than providing a quantitative error directly on the 
relative luminosity measurement, an alternative method of asymmetry calculation 
that does not rely upon measurement of the relative luminosity is used to estimate 
the uncertainty on the asymmetry values calculated using Eq. 15.11 See Section 15.6. II 
for a description of this alternative method. Refer to Appendix El for a discussion of 
potential sources of error in the determination of the relative luminosity for different 
single- and double-spin asymmetry measurements. 
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Figure 5.6: Combined asymmetry results for the two beams and particle production 
to the left and right. 

5.5.3 Fill-by-fill polarization correction 

The average beam polarization in the 2001-02 run was 15 ± 5%, with the 5% rep- 
resenting a systematic scale uncertainty, discussed below. Unpolarized protons in 
the beam act to dilute the physics asymmetries being measured. It is necessary to 
correct for this dilution, which is done by the factor of p-± — in Eq. 15.11 

-*l> cam 

The beam polarization varies fill by fill and is typically different for the two 
beams. In the 2001-02 run, the yellow beam frequently had slightly higher po- 
larization than the blue beam. The polarization values for both beams for all fills 
included in this analysis are given in Table l5~Tfl Statistical uncertainties on the beam 
polarization were on the order of 10~ 3 in absolute polarization. These uncertainties 
were negligible compared to the statistical uncertainties on the yields as well as 
compared to the systematic uncertainty on the polarization (see separate discussion 
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Fill 


-Pbluc 


ycllUW 


2222 


0.12 


0.14 


2226 


0.22 


0.23 


2233 


0.17 


0.19 


2235 


0.16 


0.24 


2244 


0.14 


0.18 


2251 


0.09 


0.16 


2266 


0.09 


0.10 


2275 


0.12 


0.14 


2277 


0.08 


0.11 


2281 


0.13 


0.11 


2289 


0.15 


0.15 


2290 


0.17 


0.21 


2301 


0.15 


0.17 


2304 


0.09 


0.16 



Table 5.6: Fill-by-fill beam polarization values for the 14 RHIC fills used, 
below) and were not incorporated into the final error on the asymmetry values. 

5.5.4 Acceptance correction 

The transverse single-spin asymmetry, An, is an azimuthal or "left-right" asymme- 
try. One can consider only particle production to the left or right of the polarized 
beam and calculate the asymmetry in production from spin-up versus spin-down 
bunches, as in Eq. 15.11 Performing the calculation in this way, detector acceptance 
effects cancel. Alternatively, one can consider only particle production from up- 
or down-polarized bunches and calculate the asymmetry in production to the left 
versus the right, as given in Eq. 15.21 

A 1 1 Nl - U ^ Nr (5.2) 

-fbeam ( | 

COS<y9|) Nl + TZaccNu 

for bunch polarization in the upward direction. Here Nl and Nr are the number 
of neutral pions produced in the left and right detector arms with respect to the 
polarized beam direction, and TZ^c = is the relative acceptance of the left and 
right detector arms. In this way, luminosity effects cancel. From Eq. 15.21 the left- 
right nature of the asymmetry is clear. However, maximal effects, i.e. the greatest 
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(\cos(p\) 


Ideal 


Actual: MB 


Actual: 2x2 


West Arm 


0.943 


0.955 


0.949 


East Arm 


0.883 


0.880 


0.874 


Both Arms 


0.909 


0.920 


0.913 



Table 5.7: Comparison of results for the average value of | cosy?| for ideal and actual 
detector acceptances. Results for the west and east arms are significantly different 
because the uppermost sector was not included in the west. Differences between the 
MB and triggered data are due to the distribution of dead or masked tiles in the 
ERT. 

and least particle production, are at 90° from the direction of the spin vector, which 
was vertical in the entire 2001-02 run. Thus integrating particle yields over the 
entire azimuthal coverage of the central arm spectrometers would lead to a dilution 
of the true physics asymmetry. The factor of y, ^r~sj in both Eq. 15. Il and l5. 21 corrects 
for this dilution. Note that ip = 0° is in the horizontal plane, implying no dilution 
at if = 0° or (f = 180°. 

The value of (| cos(p\) over the idealized azimuthal coverage of the detector arms 
can be calculated analytically by Eq. 15.31 



J \ cos <p\d(p 

<l COS(^|> = j^- (5.3) 

To account more carefully for dead areas in the EMCal and in the ERT, which 
would require detailed attention to the distribution of these areas in order to make 
an analytical calculation, the average values of | cos used in the analysis were 
calculated directly from the data using Eq. 15 A\ where j is a sum over all photon 
pairs with an invariant mass which fell under the 7r° mass peak. A comparison of 
results for (| cos (p |) determined analytically and from the data is given in Table IB~71 



(|cosy>|) 



N 

E I cosv?j| 

j=l 

N 



(5.4) 
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5.5.5 Subtraction of background asymmetry 

Subtraction of the asymmetry of the background is performed as given in Eq. 15.51 
taking into account the fraction of background under the n° mass peak. A^ ak 
indicates the asymmetry of all photon pairs falling under the 7r° peak; note that it 
has generally been written simply as A^ up until this point. Equation 15. 61 gives the 
prescription for calculation of the final statistical uncertainty on the tc° asymmetry 
after subtraction of the background asymmetry. 

At = A \ TAn (5.5) 



°" 4 peak + r2(j2 A bz 



V A nO = " (5.6) 

n 1 — r 

The same technique to handle background in a tt° asymmetry analysis was used for 
the longitudinal double-spin asymmetry [T7|. Refer back to Table IB~3l for the fraction 
of background in each pt bin. After cuts to reduce the background, it ranged from 
34% in the 1-2 GeV/c p T bin to 5% in the 4-5 GeV/c bin. 

It is not possible to measure the asymmetry of the background under the peak 
directly. Differentiation between true neutral pions and combinatorial background 
is only possible statistically and not on an event-by-event basis. Therefore it must 
be estimated in order to correct for it. The asymmetries of two different background 
invariant-mass regions were studied as estimates of the asymmetry of the background 
under the 7r° peak: 50-MeV/c 2 regions immediately around the 7r° mass peak (60- 
110 MeV/c 2 and 170-220 MeV/c 2 ), and the invariant -mass region above the 7r° but 
below the r\ (250-450 MeV/c 2 ). The final results published in ^3] subtracted the 
asymmetry of the two 50-MeV/c 2 regions. It was felt that this invariant-mass region, 
being closer to that directly under the peak, was likely to reflect the asymmetry of 
background under the peak more accurately. However, the background region used 
was found to have little effect on the final results. The similarity in the background 
asymmetries for the two different invariant-mass regions lent confidence to their 
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p T (GeV/c) 


N 


cx.bg 


1-2 


-0.007 


0.009 


2-3 


-0.031 


0.034 


3-4 


0.036 


0.123 


4-5 


0.42 


0.39 



Table 5.8: Asymmetry results of background photon pairs falling within 50-MeV/c 2 
regions around the 7r° mass peak. 



p T (GeV/c) 


N 


O" 4lr o 


1-2 


-0.005 


0.012 


2-3 


-0.012 


0.020 


3-4 


-0.016 


0.047 


4-5 


0.052 


0.109 



Table 5.9: Background-subtracted ir° asymmetries, using 50-MeV/c 2 mass regions 
around the 7r° peak as the background. 

validity in estimating the asymmetry of the background under the 7r° peak. See 
Section 15.6.51 for further comparison of the results obtained from the two different 
background regions. 

The asymmetries measured for photon pairs falling in the background invariant- 
mass region immediately surrounding the peak are shown in Table 15.81 The asym- 
metries after background-asymmetry subtraction for the two 50-MeV/c 2 regions are 
given in Table 1531 The asymmetry results before and after subtraction of the back- 
ground asymmetry are shown in Figure 15.71 It can be seen that correction for the 
background asymmetry made only a small difference in the results. The mean pt 
in each bin was adjusted to account for possible differences between the mean px of 
the true neutral pions and the background pairs under the peak. The final value is 
calculated by Eq. 15.71 

peak bg 

p T ° = Pt ~ rPT (5.7) 
1 — r 



and given for each p? bin in Table 15.101 
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Figure 5.7: Asymmetry results before and after correction for the asymmetry of 
the background in the invariant-mass regions immediately surrounding the 7r° peak. 
Background-corrected points are shifted down by 50 MeV/c from the center of the 
bin for readability. 



( PT ) (GeV/c) 


7T° peak 


Background 


Final 


1.40 


1.31 


1.45 


2.34 


2.28 


2.34 


3.35 


3.33 


3.36 


4.38 


4.37 


4.38 



Table 5.10: Mean px values of the background and of neutral pions after correction 
for background. 
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Figure 5.8: Number of photon pairs per MB event in 7r° invariant- mass range, 
1 < p T < 2 GeV/c, west arm. 

5.5.6 Calculation of statistical uncertainties 

In order to calculate the statistical error, accounting for the fact that multiple neu- 
tral pions could be produced per collision, the multiplicity distributions of 7r°'s per 
event for the different detector arms were determined. See Figures 15.81 and 15.91 
for sample multiplicity distributions in the west arm. Note that these distribu- 
tions in fact indicate the multiplicity per event of photon pairs with an invariant 
mass between 120 and 160 MeV/c 2 , thus including both real 7r°'s as well as false 
combinatorial pairs. From these multiplicity distributions, the degree to which the 
distribution is non-Poisson was calculated in a simplistic way by taking the ratio of 
the distribution's RMS to the square root of its mean. For a Poisson distribution, 
this ratio would be 1. 

Table I5~TT1 shows the degree to which the photon pair yield is non-Poisson, calcu- 
lated as <jfc/V< k >, where < k > is the mean number of pairs per triggered event 
and Ok is the RMS of this multiplicity distribution. These values were determined 
for the MB data, and then the uncertainty on each individual yield, N, was taken 
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9 10 

ji°/event 



Figure 5.9: Number of photon pairs per event in 7r° invariant-mass range, 2 < px < 
3 GeV/c, west arm. 



p T (GeV/c) 


7T° peak 


Background 




West 


East 


West 


East 


1-2 


1.08 


1.06 


1.16 


1.11 


2-3 


1.02 


1.02 


1.08 


1.05 


3-4 


1.01 


1.01 


1.04 


1.02 


4-5 


1.02 


1.00 


1.41 


1.00 



Table 5.11: Degree to which the yield is non-Poisson for MB photon pairs falling in 
the 120-160 MeV/c 2 and (60-110 or 170-220 MeV/c 2 ) invariant-mass regions. 

to be j^- \fN . The 4-5 GeV/c background bin suffers from very low statistics. 
Rather than using the calculated value of 1.41 for the west arm, 1.04, the same as 
for the 3-4 GeV/c background bin, was used. 



5.5.7 Asymmetry scale uncertainty 

In the 2001-02 run, only the pC polarimeter was available in RHIC. As stated in 
Section 14.2.41 the analyzing power in the process utilized by the pC polarimeter to 
measure the beam polarization was originally measured by AGS experiment E950 to 
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±30% • The total systematic error on the measurement of the beam polarization 
was derived from a relative systematic uncertainty on the RHIC beam measurement 
of 15%, the 30% relative uncertainty on the analyzing power of the process, and a 
relative uncertainty of 10% in the change in analyzing power from a beam energy 
of 22 GeV at the AGS to 100 GeV at RHIC [1101 ■ Adding these uncertainties in 
quadrature gave a total relative systematic uncertainty on the beam polarization 
of ±35%. This uncertainty is a scale uncertainty; it affects asymmetry values and 
statistical errors, generally proportional to 3=, in the same way, preserving the 



significance of each point from zero. 

5.6 Studies and checks 

5.6.1 Alternative asymmetry calculation 

As described above in Section I5.5.4| azimuthal transverse single-spin asymmetries 
can be considered for a single polarization direction comparing particle production 
to the left and right of the beam, or for particle production on a single side of the 
beam comparing different polarization directions. Equation l5.8l combines yields from 
up- and down-polarized bunches and from the left and right halves of the detector 
such that systematic errors are reduced. 



In particular, the acceptance and luminosity asymmetries cancel out to several or- 
ders. See |143j for a detailed discussion of this and other methods of calculation for 
transverse single-spin asymmetries. It should be noted that while Eq. 15.11 is mathe- 
matically exact, Eq. 15. 81 is an approximation, albeit an excellent one for the purposes 
of this analysis. It should also be noted that Eq. 15.81 is only suitable for transverse 
single-spin analysis, while Eq. 15.11 has an equivalent for longitudinal double-spin 
analysis (see Eq. IA.1|) . 

The asymmetry is determined for every fill, then averaged over all fills, as in the 




(5.8) 
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Pt [yjvv / c) 


1 < o q -j y i 
J_>t3CLlil 






1-2 


Blue 


-0.010 


0.012 




Yellow 


-0.001 


0.009 


2-3 


Blue 


-0.017 


0.026 




Yellow 


-0.009 


0.020 


3-4 


Blue 


-0.002 


0.068 




Yellow 


-0.032 


0.053 


4-5 


Blue 


0.06 


0.16 




Yellow 


0.08 


0.12 



Table 5.12: Asymmetry results as determined by Eq. 



calculations using Eq. 15.11 Table 15.121 shows the results obtained from Eq. 15.81 for 

triggered events. 

The results of the asymmetry calculations obtained using Eq. 15. H and Eq. 15 .81 can 
be seen together in Figure IB*. 101 The figure is for the yellow beam; results for the blue 
beam are similar. The two methods agree extremely well. The dominant systematic 
uncertainty in the results from Eq. 15.11 is expected to be from the determination 
of the relative luminosity; therefore, systematic uncertainties are calculated from a 
direct quantitative comparison of the asymmetry results obtained from these two 
methods of calculation. 



A% = A, 



A systematic uncertainty, a sys , was calculated for each bin from A l N 
a A i , and a A 2 . The index '1' refers to results from Eq. 15.11 '2' to results from Eq. 15.81 
The uncertainty on the difference was calculated following a prescription for results 
obtained by applying two different methods to data sets that are 100% correlated 
(exactly the same data) 



the correlation p 
a 



a 



sys 



In this case, a a 



IAI 



a 2 Al + a 2 



A 2 

N 



with 



a a was taken in cases where |A| > a a, and 
sys — was taken in cases where |A| < a a- See Tables 15*. 1 31 and 15 . 141 for the values 
used in the calculation and the results. The systematic error was also calculated 
as the square root of the difference of the squares, a sys = a/|A| 2 — a A , to see if it 
affected the results. At the level of 10~ 3 , which is the precision quoted in [T3], it 



did not. Final systematic uncertainties calculated using this method, averaged over 
the two beams, are shown in Table 15.151 
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Figure 5.10: Comparison of asymmetry results obtained using Eq. 15.11 (Method 1) 
and Eq. 15.81 (Method 2), shown here for the yellow beam. Points for Method 2 are 
shifted down by 50 MeV/c from the center of the bin for readability. 



Pt 
(GeV/c) 


A 1 

N 

(xlO 2 ) 


A 2 

N 

(xlO 2 ) 


A 

(xlO 2 ) 


a ai 
(xlO 2 ) 


(xlO 2 ) 


0"A 

(xlO 2 ) 


@sys 

(xlO 2 ) 


1-2 


-0.291 


-0.292 


0.001 


0.816 


0.816 


0.000 


0.001 


2-3 


-1.79 


-1.79 


0.00 


2.07 


2.07 


0.00 


0.00 


3-4 


1.61 


1.61 


0.00 


5.51 


5.51 


0.00 


0.00 


4-5 


4.23 


4.17 


0.06 


13.0 


13.0 


0.00 


0.06 



Table 5.13: Agreement of asymmetry results from the two methods of calculation, 
blue beam. Note that these values were obtained before final cuts were performed. 
See text for further explanation. 
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Pt 
(GeV/c) 


4 1 

JV 

(xlO 2 ) 


A 2 

N 

(xlO 2 ) 


A 

(xlO 2 ) 


a ai 
(xlO 2 ) 


a a2 

N 

(xlO 2 ) 


0"A 

(xlO 2 ) 


@sys 

(xlO 2 ) 


1-2 


-0.974 


-0.975 


0.001 


0.636 


0.635 


0.001 


0.000 


2-3 


-2.50 


-2.50 


0.00 


1.61 


1.61 


0.00 


0.00 


3-4 


-3.29 


-3.30 


0.01 


4.27 


4.28 


0.01 


0.00 


4-5 


6.81 


6.99 


-0.18 


10.0 


10.1 


0.1 


0.1 



Table 5.14: Agreement of asymmetry results from the two methods of calculation, 
yellow beam. The data samples are 100% correlated, so agreement at better than 
the statistical level is expected. Note that these values were obtained before final 
cuts were performed. See text for further explanation. 



Pt 
(GeV/c) 


Avg a sys 
(xlO 2 ) 


1-2 


0.000 


2-3 


0.00 


3-4 


0.00 


4-5 


0.1 



Table 5.15: Systematic uncertainty on the neutral pion asymmetry, calculated for 
each bin and averaged over both beams. 

While the comparison of the asymmetry results from the two different methods 
of calculation provided a strong check on the relative luminosity, the main identified 
potential source of systematic error present in the asymmetries as calculated by 
Eq. 15.11 a number of other checks on the results were performed. These other 
checks are described in the following sections. 

5.6.2 Left and right detector arms 

Calculating the asymmetry using Eq. 15.11 necessarily gives separate results for the 
left and right sides of the polarized beam. There is no overlap in the particle yields 
from the two detector arms; therefore, the expected agreement between the two 
results is that for uncorrelated samples. The uncertainty on the difference in results 
is taken to be <ta = , /|c 2 left + a 2 ; ht — 2pa j$.<&cr with the correlation p = 0. 

Y A N A N N N 

As can be seen from Table 15.161 the difference in the results from the two detector 
arms was within the uncertainty on the difference, and the results were in agreement. 
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p T (GeV/c) 


A left 
N 


A right 
N 


A 


a lcft 


fright 


OA 


1-2 
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0.026 
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0.0941 
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0.01 


0.13 


-0.12 


0.15 


0.16 


0.219 



Table 5.16: Evaluation of agreement of asymmetry results from particle production 
to the left and right of the polarized beam, shown for the yellow beam. A is the 
difference between the asymmetries for the two arms, a a is the uncertainty on the 
difference. 
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Figure 5.11: Comparison of results for particle production to the left and right side, 
blue beam polarized. Points for the right of the beam are shifted down by 50 MeV/c 
from the center of the bin for readability. Statistical agreement is expected. 



Results for both detector arms can be seen together for the blue and yellow beams 
in Figures 15.111 and 15.121 respectively. 
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Figure 5.12: Comparison of results for particle production to the left and right side, 
yellow beam polarized. Points for the right of the beam are shifted down by 50 
MeV/c from the center of the bin for readability. Statistical agreement is expected. 

5.6.3 Two independent beams 

While the results for the blue and yellow beams use the same events and yields, 
they are combined in a different way, taking into account the spin direction of 
either one beam or the other, leading to (nearly) statistically independent mea- 
surements. As the correlation between results from the two beams is believed to 
be small but is unknown, it is assumed to be zero in evaluating the agreement of 
the results. Thus as for the case of comparing results from the two detector arms, 
the uncertainty on the difference of the results from the two beams is taken to be 
°a = w 1 0" ^ blue + °" 2 ,yciiow — 2pa A biuc(j A ydiow |, with p = 0. It can be seen in Table I5~TTI 

Y A N N N N 

that the results from the two beams agree as expected. Results for both beams can 
be seen together and evaluated by eye in Figure 15.131 
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Table 5.17: Evaluation of agreement of asymmetry results from blue and yellow 
beams. A is the difference between the asymmetries for the two beams. <ta is the 
uncertainty on the difference. 
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Figure 5.13: Comparison of results obtained for the two polarized beams. Points 
for the yellow beam are shifted down by 50 MeV/c from the center of the bin for 
readability. Statistical agreement or better is expected. 
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a A N 
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bg2 


bgl 


bg2 


1-2 


-0.008 


-0.014 


0.005 


0.004 
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-0.006 


0.008 
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0.013 


3-4 


-0.012 


0.015 


0.079 


0.055 


4-5 


0.00 


0.03 


0.21 


0.14 



Table 5.18: Asymmetry results of photon pairs falling within the two 50-MeV/c 2 
regions around the mass peak (bg 1) and within 250-450 MeV/c 2 (bg 2). Note that 
this check was performed before the final cuts on the data sample; thus, the results 
from bg 1 shown here differ from the final background results. 

5.6.4 Triggered and minimum-bias data 

Results for the physics asymmetry, An, can also be compared for the triggered and 
minimum-bias data samples. In this case the correlation between the samples is 
poorly understood. In the higher p? bins, where the trigger was more efficient, 
nearly all pions in the MB sample should have fired the trigger and been present 
in the triggered sample as well, making the MB sample nearly a direct subset of 
the triggered sample. In the lower bins, this correlation should be lower but still 
significant. Because of the unknown correlation and the fact that the MB sample 
is severely statistically inferior to the triggered sample, no direct evaluation of their 
agreement was performed. However, it can be seen in Figure 15*. 141 that the MB and 
triggered results do not exhibit notable disagreement. 

5.6.5 Different background regions in invariant mass 

As described above, the asymmetry of two different background mass regions in 
invariant mass was investigated: 50-MeV/c 2 regions around the 7r° mass peak (60- 
110 MeV/c 2 and 170-220 MeV/c 2 ) and the mass region between the 7r° and rj (250- 
450 MeV/c 2 ). Table [5~TH1 shows the background asymmetries as calculated by Eq. 15. II 
for the two background regions. The asymmetries are similar. The background 
asymmetries are consistent with zero for both background regions for 2 < pt < 5. 
In the 1-2 GeV/c pt bin, both regions suggest a slightly negative asymmetry. 

In Figure IK. 151 a direct comparison of the asymmetry results after subtraction of 
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Figure 5.14: Comparison of results for minimum-bias and triggered events, shown 
here for the yellow beam and east detector arm. Points for the triggered data are 
shifted down by 50 MeV/ c from the center of the bin for readability. The correlation 
between the samples is unknown. 

the two background asymmetries is shown. The background region used has little 
effect on the final asymmetry. 

5.6.6 Different neutral pion invariant-mass integration re- 
gions 

As an additional check on the sensitivity of A^ to the background under the 7r° peak, 
three different integration regions for the it mass were examined: 120-160 MeV/c 2 , 
110-170 MeV/c 2 , and 100-180 MeV/c 2 . Figure IB. 161 shows a comparison of asym- 
metries obtained for peak integration from 120-160 MeV/c 2 and 100-180 MeV/c 2 . 
There is little effect on the result from the amount of background included, pro- 
viding additional evidence that the background under the peak does not affect the 
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Figure 5.15: Comparison of asymmetries obtained after subtracting the asymmetry 
of two different background invariant-mass regions, "bg 1" indicates the asymmetry 
after subtraction of the asymmetry around the peak; "bg 2" is after subtraction 
of the asymmetry of the 250-450 MeV/c 2 background region. Points for bg 2 are 
shifted down by 50 MeV/c from the center of the bin for readability. Note that this 
check was performed before the final cuts on the data sample; thus the results differ 
from the final ones. 

asymmetry greatly. 



5.6.7 Bunch shuffling 

A technique called " bunch shuffling" can be utilized to check for uncorrelated bunch- 
to-bunch and fill-to-fill systematic errors. For each bunch crossing the spin direction 
is reassigned randomly, and then the new asymmetry with false spin dependence is 
recalculated. This procedure is repeated many times. With random reassignment 
of the spin direction to each crossing, a Gaussian asymmetry distribution centered 
around zero is expected. One must take care in choosing the exact procedure used for 
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Figure 5.16: Comparison of results for different 7T° peak integration regions. Points 
for the wider integration region are shifted down by 50 MeV/c from the center of 
the bin for readability. 



this study. If done correctly, the root-mean-square (RMS) width of the "shuffled" 
(non-physics) asymmetry distribution should correspond to the statistical uncer- 
tainty on the physics asymmetry, since fluctuations in the calculated asymmetry 
should not be due to any spin dependence but rather to statistical fluctuations in 
event-by-event particle production. If the RMS width of the shuffled asymmetry 
distribution is larger than the statistical uncertainty on the physics asymmetry, it 
should reflect the presence of elements that broaden the distribution beyond statis- 
tical fluctuations in particle production, e.g. some source of bunch-to-bunch system- 
atic error in the physics asymmetry. It should be noted, however, that the expected 
quantitative agreement between the statistical uncertainty and RMS width of the 
shuffled distribution is not entirely understood. As is discussed in Section I5.fi.7.1| 
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a study done for this thesis found that it is also possible to obtain widths smaller 
than the statistical uncertainties. 

It is easier to consider the validity of various procedures assuming large asym- 
metries, for example 1 (100%) . Taking one bunch at a time and randomly assigning 
its spin direction to be up or down, with no further constraints, it is theoretically 
possible to make the spin assignments exactly as the true, original spin directions 
(up phys -> up shuf , down phys -> down shuf ), or exactly opposite (up phys -> down shuf , 
down p h ys — > up shuf ), yielding a shuffled asymmetry distribution that can range from 
-1 to +1, even with a wealth of statistics. This method would generally give dis- 
tribution widths wider than the statistical uncertainties on the physics asymmetry 
values. Selecting any half of the bunches to be assigned up and the other half to be 
assigned down creates a similar situation for true relative luminosity values close to 
1, i.e. a nearly equal number of events coming from bunches with spin up and spin 
down: it is possible to assign the shuffled spins to be nearly the same as the original 
physics spin directions. 

The procedure utilized in this analysis was to assign half of the original up spins 
to down and keep the others as up, and similarly assign half of the original down spins 
to up while keeping the others down. Thus the particle production gets redistributed 
evenly, at the crossing (not event) level, between up and down spins, and even 
for large physics asymmetries, shuffled asymmetry will always be (nearly) zero. 
Repeating this procedure many times should yield a distribution around zero that 
is due to event-by-event statistical fluctuations in particle production (i.e. 0, 1, 2, 
etc. neutral pions produced in a given event, with a certain probability distribution). 

For this analysis, Eq. 15.81 was used to calculate the shuffled asymmetries, as it 
avoided the complication of recalculating the relative luminosity for each iteration. 
The bunch shuffling procedure was performed 1000 times. Refer to Figure 15.171 
for examples of shuffled asymmetry distributions from the data set used for this 
analysis; the expected increase of the distribution width for the smaller-statistics 
(higher-pr) bins can be clearly seen. In Table l5~T9l it can be seen that the mean of 
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Figure 5.17: Bunch- shuffled asymmetry distributions for triggered data, shown here 
for the blue beam. 

the shuffled asymmetry distributions was zero at the level of 10~ 4 or better. The x 2 
distributions for a fit to a constant asymmetry across all shuffled fills can be seen 
in Figure 15.181 The dotted line shows the expected \ 2 distribution for 13 degrees 
of freedom (asymmetry values for 14 fills, fit to a constant), which are in agreement 
with the data. 

Table 15.201 shows the statistical uncertainties on the physics asymmetries com- 
pared to the RMS widths of the shuffled asymmetry distributions for the blue and 
yellow beams. In order for the uncertainties on the physics asymmetries to reflect 
simply statistical distributions in particle production, the beam polarization was not 
factored in. For this reason, the statistical uncertainties are exactly the same for 
the blue and yellow beams, the results for each using the same total set of events. 
The % difference is calculated as 100 
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Figure 5.18: Bunch-shuffled \ 2 distributions from triggered data, shown here for the 
blue beam. The dashed lines indicate the expected distributions. 
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Table 5.19: Mean values of shuffled asymmetry distributions. 
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Table 5.20: Comparison of statistical uncertainties on the physics asymmetries 
(o^phys) and RMS widths of shuffled asymmetry distributions (a A ^t). 

The shuffled widths are systematically smaller than the statistical uncertainties 
calculated for the physics asymmetry. This outcome was at first a surprise, as the 
expectation was that the RMS width should have been the same as or greater than 
the physics statistical uncertainty. However, a subsequent Monte-Carlo study of 
the bunch shuffling technique, discussed below, corroborated the tendency for the 
width of shuffled distributions for a zero physics asymmetry to be narrower than 
the statistical uncertainty. 

5.6.7.1 Bunch shuffling Monte Carlo 

A Monte Carlo study of the bunch shuffling technique was performed to better quan- 
titatively understand the expected agreement between the RMS widths of shuffled 
asymmetry distributions and the statistical uncertainties of the physics asymmetry 
values, for different simulated physics asymmetry values. 

In the Monte Carlo study, 14 fills were assumed, the same number as in the 
actual data set. In order to simulate the variation in the total particle yield per 
fill, the actual total tc° yields per fill for the 1-2 GeV/c px bin were used as the 
simulated yields per fill. The number of bunches per beam was taken to be 48, 
selected because it was very close to the typical value of 46 in the actual data set, 
and because it was divisible by four, so that it was possible to assume half of the 
bunches were spin up and half were spin down, then easily reassign exactly half in 
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each spin group to the wrong spin direction. The desired physics asymmetry to 
simulate was selected, which determined the average yield per bunch with spins up 
and down in each fill. Simulated yields for each bunch crossing were produced by 
sampling from Poisson distributions around these asymmetry-dependent averages. 
In this way, the uncertainty on the yields themselves was purely Poisson and well 
understood. The uncertainty on the simulated asymmetry was purely statistical 
and calculated by performing straightforward error propagation on the asymmetry 
formula assuming only uncertainties on the yields; there were no effects due to 
detectors, triggers, polarization measurements, or any other factors incorporated. 
The simulated yields were then input into the same software program written to 
handle the bunch shuffling of the data. The RMS width of the resulting shuffled 
asymmetry distribution was then compared to the statistical error on the simulated 
"physics" asymmetry. 

The effects of varying the simulated physics asymmetry, the total statistics, 
and the number of bunches per beam were studied. In Table 15.211 the results for 
two different statistical sample sizes as well as three different simulated asymmetry 
values are shown. There are several items to note here. Despite incorrectly assigning 
the spin direction for exactly one half of the bunches, the width of the shuffled 
distribution increases noticeably as a function of the simulated physics asymmetry. 
The statistical uncertainty on the physics asymmetry is itself weakly but inversely 
dependent on the physics asymmetry, leading to a shuffled width that is more than 
double the simulated uncertainty for a 90% simulated asymmetry. For the case of a 
zero simulated asymmetry, corresponding to the results of the present analysis (see 
Section I5~7j) . a shuffled width smaller than the simulated uncertainty is obtained for 
a statistical sample equivalent to the 7r° yield in the 1-2 GeV/c pt bin in the present 
analysis. For ten times this statistical sample, the study suggests that the shuffled 
widths tend to increase with respect to the smaller sample, for at least the zero and 
30% simulated asymmetry values. In order to investigate the possibility that the 
smaller width for zero simulated asymmetry was due to correlations in obtaining 
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shuffled asymmetries many times for a limited number of bunch crossings, a study 
with only 10 bunches per beam was performed, the results of which are given in 
Table E221 Performing 10,000 reassignments of the spin direction of five out of ten 
bunches should have led to many duplicate reassignments and thus duplicate shuffled 
asymmetry values in the distribution. (Note that the difference in the number of 
iterations in the data versus the simulation was due simply to technical reasons.) 
For the smaller statistical sample size, the shuffled distribution width was only 
slightly smaller for 10 bunches per beam than 48 bunches; for the larger sample 
size the reduction was more significant. While the degree of correlation should 
not depend on the value of the physics asymmetry, it may be that the reduction 
of the shuffled distribution width due to correlations only makes it less than the 
statistical uncertainty in the case of zero asymmetry because for larger asymmetries, 
the dependence of the width on the asymmetry is a larger effect. 

As an additional means of investigating correlations as the cause of shuffled 
distribution widths being smaller than the statistical uncertainties, an attempt was 
made to remove all correlations. To achieve this, the yield for each crossing was 
resampled for each iteration of the reassignment of spin directions. This procedure 
effectively incorrectly assigned the spin direction for the yields from the bunches 
in many different fills (Jj effective fills = (j shuffling iterations), all having the same 
physics asymmetry, rather than repeatedly reassigning the spin direction for the 
yields in a single fill or small number of fills. The results of this exercise are shown 
in Table I5.23| the shuffled width is within 1% of the statistical uncertainty on 
the simulated physics asymmetry, suggesting that the smaller shuffled distribution 
widths may be due to correlations. Note that this procedure is not applicable to 
real data, due to the fact that there is always a limited number of fills. 

While the Monte Carlo investigation confirmed the pattern seen in applying 
the bunch shuffling technique to the actual data, the merits and limitations of the 
technique are still not completely understood. Improved procedures for randomly 
reassigning the spin direction of the bunches in order to obtain better agreement 
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Table 5.21: Bunch shuffling results for two different statistical sample sizes and 
three different simulated asymmetry values, assuming 48 bunches per beam. The 
"Statistics" column indicates the factor times the 1-2 GeV/c px bin 7r° yields in the 
actual data sample. 10,000 shuffles were performed. 
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Table 5.22: Comparison of simulated shuffling results for 10 and 48 bunches per 
beam. 10,000 shuffles were performed. 

between statistical uncertainties on the physics asymmetries and shuffled asymmetry 
distribution widths may exist. A more thorough future study would be valuable. 
This method of checking for uncorrelated bunch-to-bunch and fill-to-fill systematic 
errors is particularly important for double-spin asymmetry measurements, in which 
many of the checks available to a single-spin analysis are not available. 

5.7 Results 

Final asymmetry results for mid-rapidity neutral pions from 200-GeV polarized p+p 
collisions for 1 < p T < 5 GeV/c, as published in [15] . are given in Figure l5~T9l and 
Table IS~2"4l They utilized the triggered data sample and the 50-MeV/c 2 regions in 
invariant mass around the 7r° peak for the background correction. The asymmetries 
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Table 5.23: Modified bunch shuffling simulation with correlations removed (see text). 
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Figure 5.19: Final mid-rapidity neutral pion transverse single-spin asymmetry. The 
error bars represent statistical uncertainties. 
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Table 5.24: Neutral pion transverse single-spin asymmetry values and statistical 
uncertainties for all photon pairs falling within the ir° mass peak, for the background 
(bg), and for the tt° background-corrected. The third column (r) indicates the 
background contribution under the 7r° peak. An A^ scale uncertainty of ±35% is 
not included. 

are consistent with zero within a few percent for all px bins. For further discussion 
of these results and their implications, see Chapter UJ 
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0.5-1 
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Table 5.25: Charged hadron transverse single-spin asymmetry values and statistical 
uncertainties. An A N scale uncertainty of ±35% is not included. 

5.8 Comparison to charged hadron asymmetry re- 
sults 

A similar analysis of the transverse single-spin asymmetry, A^, of inclusive charged 
hadrons at mid-rapidity for 0.5 < pt < 5.0 GeV/c was performed by F. Bauer for 
the 2001-02 data. The results of the neutral pion and charged hadron measure- 
ments have been published together in [T5] . For more details on the charged hadron 
analysis, refer to the publication. Charged hadron asymmetry results were obtained 
separately for positively and negatively charged particles; they are presented in 
Table IB~2*5l While a clear charge dependence has been observed in transverse single- 
spin asymmetries for forward production (refer back to Section 12 .2 .5|) . the results for 
mid-rapidity neutral pions and both charges of hadron were all found to be similarly 
consistent with zero. Figure 13.201 shows the observed asymmetries for both neutral 
pions and charged hadrons together. 
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Figure 5.20: Mid-rapidity neutral pion and charged hadron transverse single-spin 
asymmetries versus mean p?. Points for positive hadrons have been shifted down by 
50 MeV/c to improve readability. The error bars represent statistical uncertainties. 
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Chapter 6 

Future prospects for transverse 
spin physics 

The experiments at RHIC have the potential to perform a wide variety of studies 
related to the transverse spin structure of the proton. Several measurements be- 
yond those completed have already been planned, and others are being proposed 
and considered. Because there are currently a number of possible mechanisms un- 
der discussion that may contribute to the large transverse SSA's observed, seeking 
measurements sensitive to or dominated by a single mechanism is especially valu- 
able. In this chapter, potential future measurements at RHIC are discussed, with 
particular focus on possible measurements by the PHENIX experiment. 

6.1 BRAHMS 

The BRAHMS experiment does not have spin rotator magnets around its IP and 
therefore does not have access to longitudinally polarized collisions. Consequently, 
its polarized proton program is completely focused on transverse spin physics. BRAHMS 
was designed to specialize in the measurement of identified charged particles at for- 
ward and mid-rapidity. It is in a unique position at RHIC to make a precision 
measurement of for identified charged particles (7r + /7r~, K + /K~, p/p) in the 
forward and backward directions with respect to the polarized beam. Preliminary 
results from pions and protons are already available from the 2004 data |14fij ; the 
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first kaon asymmetries as well as more accurate pion and proton asymmetries are 
expected from the data obtained in 2005. The kaon asymmetries and their charge 
dependence will provide additional information directly in support of or against the 
hypothesis that the large asymmetries observed are due to valence quarks. Assum- 
ing they are, positive kaons, with a quark content of us, would be expected to have 
a positive asymmetry, and negative kaons, comprised of us, would be expected to 
have an asymmetry of approximately zero or only slightly negative. More precise An 
measurements will also help to distinguish among various model calculations that 
are consistent with the present results and to put an upper bound on the Sivers dis- 
tribution function. Now that the BELLE experiment has made a first measurement 
of the Collins FF pQ, transverse SSA's will also be able to constrain transversity. 
The 2005 BRAHMS results for pions are expected to cover 0.15 < \xp\ < 0.35; lower 
beam energies in future running would allow even higher xp values to be reached. 
Current results from STAR suggest that the 7r° asymmetry may have a maximum at 
xf ~ 0.5 or 0.6 and decrease quickly beyond that; see below. Lower-energy running 
would permit exploration of the behavior of the charged-particle asymmetries as a 
function of Xp and provide information on transversity at different x values. RHIC 
is capable of colliding polarized proton beams at center-of-mass energies as low as 
50 GeV. 

6.2 STAR 

The STAR detector has the largest acceptance among the RHIC experiments. Its 
main barrel, a time-projection chamber for tracking charged particles surrounded 
by an electromagnetic calorimeter, has full azimuthal coverage for \q\ < 1, making 
it particularly suitable for multi-particle correlation measurements and jet recon- 
struction. In 2003 Boer and Vogelsang proposed a single transverse-spin di-jet mea- 
surement that could isolate the Sivers effect and probe the gluon Sivers function 
[6*H] . A non-zero Sivers function implies a spin dependence in the distributions 
of the partons within the proton, which would lead to an observable spin-dependent 
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asymmetry in A(p of (nearly) back-to-back jets. A discussion of how such an anal- 
ysis would be performed by STAR is given in [^Jj; there will likely be a significant 
result after the 2006 run at RHIC. 

As a follow-up to earlier neutral pion An measurements with a forward 7r° detec- 
tor covering 3.3 < 77 < 4.2 f3\ 1126j . STAR intends to make a precision measurement 
of the 7T° transverse SSA over a large xp range. The current results suggest that 
the asymmetry may decrease in magnitude for xf > 0.6; additional data should 
make this trend clearer, if it exists. As is the case for BRAHMS measurements, 
further A^ results from STAR will be able to constrain the Sivers and transversity 
distributions. 

Through recent upgrades, STAR now has electromagnetic calorimetry covering 
— 1 < rj < 2, allowing observation of neutral pions and direct photons within this 
acceptance. Measurement of A^ for direct photons has been proposed as a probe 
of the Sivers function for gluons |139j . As mentioned previously, direct photon 
production is dominated by quark-gluon Compton scattering (q + g — > 7 + X) 
over a wide range in photon transverse momentum at RHIC. The fact that there is 
no fragmentation involved in direct photon production eliminates contributions to 
the asymmetry from the Collins effect. Alternatively or additionally, STAR could 
measure correlated photon-jet pairs, which would offer better understanding of the 
kinematics of the direct photon production, i.e. the x values of the two partons. Jet 
production, similar to direct photon production, does not consider fragmentation to 
particular hadrons and is insensitive to the Collins effect. SSA's of photons and jets 
in events with correlated photon-jet pairs would thus access the gluon and quark 
Sivers functions, respectively, with some ability to identify the x values at which 
these functions are probed. 

6.3 PHENIX 

As discussed in Section |4~3*1 PHENIX specializes in the measurement of leptonic and 
photonic probes and high-py particles over a limited acceptance. From the modest 
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transverse-spin data sample taken in 2005 (0.16 ptT 1 , ~ 48% average polarization), 
PHENIX has already begun analysis to obtain improved mid-rapidity An results for 
neutral pions and charged hadrons, which are expected to provide tighter constraints 
on the gluon Sivers function. Future higher-statistics data samples for these particles 
at mid-rapidity will provide greater sensitivity to transversity and the Collins effect. 
See Chapter [7| for further relevant discussion. 

There is also analysis underway to obtain first results for An of single muons, 
largely from open charm decay but with significant contributions from light-hadron 
decays. The current xp reach for this measurement is ~ 0.15, lower than the xp 
values at which significant asymmetries have been observed for other particle species; 
higher xp values would become accessible with lower-energy running. A forward 
hadron An measurement using the PHENIX muon spectrometers (1.2 < \r]\ < 2.4) 
may also be possible using muons from hadron decays and the charged hadrons 
themselves that punch through the steel absorber in front of the muon tracker; 
however, careful studies will need to be done in order to understand the particle 
ratios in this sample. 

In 2006, PHENIX intends to perform a transverse single -spin di-hadron measure- 
ment similar to the STAR di-jet measurement isolating the Sivers effect, following 
the proposal in 68J. This analysis would study the spin dependence of the azimuthal 
angle between nearly back-to-back 7r°-hadron pairs, triggering on a decay photon 
from the tt° in order to obtain a higher-statistics sample. Although dilution of the 
effect is anticipated for hadron rather than jet pairs, studies have shown that the ef- 
fect should still be measurable. In examining back-to-back hadrons rather than jets, 
fragmentation to the final-state hadrons must be considered, and some sensitivity to 
the Collins mechanism is introduced. However, as described in Chapter there is 
a large contribution from gluon fragmentation to tt° production for pp < 5 GeV/c, 
to which the Collins mechanism does not apply. The effect is expected to be maxi- 
mized when the initial spin of the proton is in the same direction as the final-state 
jet axis; thus, to match the PHENIX central arm acceptance, data will be taken 
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with radially polarized collisions. 

Also similar to possible measurements at STAR, PHENIX could measure An 
of direct photons, sensitive to the Sivers function. Future upgrades extending the 
azimuthal coverage for tracking to 27r in the inner region and adding forward elec- 
tromagnetic calorimetry (0.9 < < 3.0) are expected to expand the coverage for 
this measurement as well as make 7-jet and jet-jet measurements feasible. 

It has been proposed to study the transverse SSA of D meson production at 
RHIC as a measurement of the gluon Sivers function (23 EH] • This measurement 
would be significant at any point over a wide range of xp values (—0.2 < xp < 0.6). 
D mesons are produced principally via the reaction g + g —> c+c, with contributions 
from q + q — > c + c becoming important only at very large Xf (xp > 0.6). In 
neither process is the final core polarized, excluding the Collins mechanism from 
contributing to any asymmetry that may be observed. Any transverse SSA seen for 
mid- to moderate rapidity D production would thus be a direct indication of a non- 
zero gluon Sivers function. PHENIX is currently capable of measuring open charm 
decays statistically via single electrons in the central arms and single muons in the 
muon spectrometers. In the future, a silicon vertex detector upgrade will make 
it possible to identify D mesons event-by-event. A silicon-pixel and silicon-strip 
barrel detector, covering the central arm pseudorapidity region and 2tc in azimuth, 
is expected to be installed in 2009. A silicon-strip endcap detector, covering the 
muon arm acceptance, has been proposed for 2011. 

The flavor separation of the Sivers function for u, d, u, and d quarks via An 
of forward or backward W boson production has been suggested |1371| . PHENIX 
already has a W physics program planned for the future, once 500-GeV polarized 
collisions are achieved by RHIC, making W measurements possible. The processes 
of interest will be: 

u + d -> W + -> /i + + 
d + u — > W~ — >• n~ + Up 

PHENIX will observe the final-state muons; an upgrade to trigger on the highest-pr 
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muons, which will come principally from W decays, is expected to be installed in 
2009. The muon trigger upgrade will also improve the pattern recognition in the 
muon arms, reducing background and making measurements such as single muons 
from D decays, or muon pairs from resonance decays or the Drell-Yan process (q + 
q — > £ + + £~), cleaner. 

T T T I 

The double transverse-spin asymmetry, A TT = ^ T T+gU , is another observable sen- 
sitive to transverse spin quantities. Att for the Drell-Yan process provides direct 
access to the transversity distribution. The transversity distributions for the quark 
and the antiquark provide the necessary convolution of two chiral-odd functions to 
be an allowed process in QCD, with no fragmentation or final-state interactions 
involved. Although this asymmetry is expected to be at the sub-percent level for 
\fs = 200 GeV, it could reach several percent for ^/s < 100 GeV. PHENIX already 
has an effective di-muon trigger that would be suitable for measuring Drell-Yan pairs; 
however, studies would need to be done to understand the current backgrounds. The 
muon trigger upgrade mentioned above would provide a cleaner sample. Investiga- 
tion would be necessary to optimize the beam energy in order to obtain the best 
compromise between luminosity and the size of the predicted asymmetry. A first 
direct measurement of transversity would be an exciting milestone in the field of 
transverse nucleon spin structure. 
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Chapter 7 
Conclusions 

The transverse single-spin asymmetry, An, for neutral pion production at xf ~0.0 
for 1 < px < 5 GeV/c from polarized proton-proton interactions at ^/s = 200 GeV 
has been presented. It is consistent with zero within statistical errors of a few per- 
cent. The measurement, together with a similar measurement for charged hadrons, 
represents the first mid-rapidity An result at high and collider energies; both 
have been published in Physical Review Letters JT^]. The measurement is in a 
kinematic regime where the theoretical framework of perturbative QCD has been 
demonstrated to describe the polarization-averaged cross sections for neutral pions 
and charged hadrons well for pt > 1.5 GeV/c. A pQCD framework is expected to be 
applicable in the interpretation of the polarized results down to similar transverse 
momentum values. 

The present result is consistent with the mid-rapidity results for neutral pions 
from p + p collisions at ^/s = 19.4 GeV |Hj. The measurement is complementary to 
the measurement of An for forward neutral pion production in p + p collisions at 
y^s = 200 GeV [9J, which observed asymmetries reaching ~ 30%. 

Neutral pion production at forward rapidity is expected to originate from par- 
tonic processes involving valence quarks (x > 0.1), whereas the particle production 
at mid-rapidity presented here is dominated by gluon-gluon and quark-gluon pro- 
cesses (x < 0.1). Figure mi shows the fractional contribution of different partonic 
scattering processes to the production of neutral pions at mid-rapidity at 200 GeV, 
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Figure 7.1: Relative fractional contributions of partonic processes to mid-rapidity 
7T° production at y/s = 200 GeV, calculated by W. Vogelsang. 

calculated by W. Vogelsang. As evident from the figure, 7r° production in the trans- 
verse momentum range covered by the measurement presented here is approximately 
half from gluon-gluon scattering and half from gluon-quark scattering. As such, the 
asymmetry is not very sensitive to mechanisms involving quarks, e.g. the Collins 
effect. 

Independently of the suppression or dilution of the Collins effect in the present 
results due to a dominance of gluon scattering, it has been stated by Anselmino et al. 
that the Collins effect is suppressed more generally [Ho*! . They argue that when 
all partonic intrinsic motion is taken into account in the formalism without sim- 
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plifications or approximations, the complicated azimuthal angle dependencies and 
numerous different phases involved in the Collins mechanism lead to strong sup- 
pression of the final asymmetry. They demonstrate that with saturated bounds on 
the non-perturbative components (the transversity distribution function and Collins 
FF), the Collins mechanism alone is insufficient to explain the large asymmetries 
observed in forward charged pion production at y/s = 19.4 GeV |6j. There has been 
some criticism of the formalism on which their argument is based (see for example 
the alternative treatment presented in jH]), but a consensus in support of its va- 
lidity or repudiating it has yet to be reached within the theoretical community. In 
the framework of Anselmino et at, the PHENIX mid-rapidity An results, i.e. the 7r° 
measurement presented here as well as the charged hadron measurement discussed 
in Section I5"%1 provide strong constraints on the Sivers distribution function for glu- 
ons and indicate that it is small. A publication providing a quantitative constraint 
on the gluon Sivers distribution from the PHENIX data is forthcoming [HI] . 

Regardless of more detailed interpretations that may become available, the re- 
sults are consistent with the naive pQCD expectation that transverse single-spin 
asymmetries are suppressed at high and mid-rapidity |112| \1'A'2\ . 

The transverse SSA results from the E704 experiment at Fermilab [BJ IHj and 
from the RHIC experiments jUJ 1146} ITS] are strikingly similar despite an order of 
magnitude difference in center-of-mass energy. In the forward direction (xf > 0), 
significant asymmetries were observed at both energies for the production of charged 
and neutral pions. In both similar, clear dependence of the asymmetry on the 

particle charge was observed: positive for positive and neutral pions and negative for 
negative pions, with m \A N |. Mid-rapidity (xp ~ 0) neutral pion results are 

consistent with zero at both energies. Asymmetries consistent with zero have been 
measured not only at mid-rapidity but also in the backward direction (xf < 0) at 
RHIC. The results for forward, mid-rapidity, and backward particle production 
considered all together suggest that the large SSA's observed in the forward region 
are due to valence quarks. 
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Interest in transverse nucleon-spin structure and transverse SSA's rose sharply 
when the first large asymmetries were observed in the late 1970's, contradicting 
naive expectations from pQCD. Despite similar effects observed by more than one 
subsequent experiment in the 1980 's, it was not until higher-energy data became 
available starting in the 1990's that a pQCD framework became potentially applica- 
ble in interpreting the effects, and the theoretical ideas most widely accepted today 
began to develop. 

The study of transverse nucleon-spin structure has progressed rapidly both the- 
oretically and experimentally over the course of the last several years. Notable 
experimental contributions have come from the HERMES and COMPASS polarized 
deep-inelastic scattering experiments, the BELLE experiment studying e + +e~ anni- 
hilation, and the STAR, PHENIX, and BRAHMS collaborations studying polarized 
p + p collisions at RHIC. 

The recent measurement of the Collins FF for pions by BELLE Ij represents a 
turning point in the study of the transverse spin structure of the proton. The non- 
zero result means that the Collins mechanism remains as a possible origin of the 
large transverse single-spin asymmetries observed. Moreover, it allows access to the 
transversity distribution function through single-spin asymmetries for the first time. 
Non-zero asymmetries already measured by DIS and hadronic-collision experiments 
can now be revisited and reinterpreted to obtain first constraints on transversity, 
assuming the asymmetries contain contributions from the Collins effect. 

The startup of RHIC as the world's first polarized-proton collider in late 2001 
ushered in a new era in the study of nucleon spin structure. RHIC holds great poten- 
tial for in-depth exploration of both the transverse and longitudinal spin structure 
of the proton. To date, two transverse spin publications |HEI] and one longitudinal- 
spin publication have come out of the major experiments at RHIC. Additional 
preliminary results are also available, as summarized in Section 12.21 

The transverse single-spin asymmetry for neutral pions presented here represents 
an early measurement in a rigorous program to study transverse proton spin struc- 
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ture at hard scales using a pQCD framework at RHIC. Conclusively explaining the 
large transverse single-spin asymmetries, which have been observed over an exten- 
sive range of energies, would help to link proton structure at soft and hard scales. 
A number of planned and proposed future transverse spin measurements have been 
described in ChapterEl Some of the measurements proposed would be able to isolate 
particular mechanisms. These measurements will be essential in order to disentangle 
the several possible contributions to the large observed asymmetries currently under 
discussion and to begin to understand the transverse spin structure of the proton. 
It should be possible to measure transversity, the last remaining leading- twist, k?- 
integrated distribution function that is completely unknown, and to start to measure 
the various transverse-momentum-dependent (TMD) distribution functions such as 
the Sivers function. 

As discussed in Section 12.2.51 TMD distribution functions are related to the 
orbital angular momentum of partons within the proton. Thus measurement of 
the Sivers function could shed light on parton OAM, which remains to date nearly 
inscrutable. No proposal for a clear and direct experimental measurement of OAM 
has yet been set forth. Therefore, measurement of the Sivers function could provide a 
starting point for elucidating this still-opaque aspect of nucleon angular-momentum 
structure. 

Despite the fact that the proton is one of the most commonplace and stable 
particles in existence, a fundamental component of ordinary matter, the path to 
unraveling proton structure has been a long one, and the journey is not yet finished. 
The final picture will not be complete without full description of its momentum, 
helicity, and transverse-spin structure. A comprehensive understanding of the pro- 
ton, in many ways the embodiment of QCD, implies an understanding of the strong 
force, one of the four fundamental forces in nature and the foundation for all of 
nuclear physics. 
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Appendix A 

Relative luminosity considerations 

In a collider environment, different considerations must be made in determining 
the relative luminosity of crossings with different spin configurations in the case of 
double- or single-spin asymmetry measurements, and for transverse or longitudinal 
spin in the case of single-spin asymmetry measurements. Care must be taken in 
designing a MB trigger that has no spin-dependent bias and can make an accurate 
measurement of the relative luminosity. To make double-spin asymmetry measure- 
ments, the principal concern is that there may be a non-zero double-spin asymmetry 
in the production of particles that fire the MB trigger. For PHENIX, the MB trigger 
fires when at least one charged particle produces a hit in each BBC (3.0 < 77 < 3.9). 
A physics asymmetry in the production of such charged particles in the kinematic 
range to fire the BBC's would thus lead to a spin-dependent bias in the "minimum- 
bias" trigger. While it is extremely difficult to prove that there is no such bias, 
checks can be performed against alternative "minimum-bias" triggers, sensitive to 
different physics processes. The relative luminosity measured by the different pro- 
cesses can be compared. This technique has been used in the measurement of the 
double-longitudinal asymmetry in neutral pion production at PHENIX |17j . 

A MB detector used simultaneously as a relative luminosity detector will bias the 
measurement if there is a double-longitudinal asymmetry in the particles to which it 
is sensitive. If the detector is capable of measuring the multiplicity of the produced 
particles, it may be possible to understand this bias. This will not be the case if it 
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APPENDIX A. RELATIVE LUMINOSITY CONSIDERATIONS 



has only a simple binary hit/no-hit response. 

The physics asymmetry in the MB/relative luminosity detector could also poten- 
tially be measured by looking for dependence of the relative luminosity on beam po- 
larization. If for example the same-helicity cross section is higher than the opposite- 
helicity cross section in the case of double-longitudinal observables, then the mea- 
sured relative luminosity should be enhanced when a higher fraction of the beam is 
polarized. 

An additional idea to measure the physics asymmetry in the MB/relative lumi- 
nosity detector would be to examine the MB rate change upon changing the spin 
combinations in a single fill. Start for example with all bunches in both beams hav- 
ing positive helicity (designated by "+"), creating only ++ collisions, and measure 
the MB event rate. Use the spin flipper described in Section r4.2.fil to flip all spins 
in a single beam in order to achieve H — collisions, and measure the MB event rate 
again. Assuming the spin flipper does not affect other beam conditions, a different 
MB rate implies a different physics cross section for ++ and H — helicity combina- 
tions. This procedure could be repeated several times throughout a fill. In order to 
check for other beam effects, it would be possible to compare +- 1- event rates after 
zero and two beam flips, or to compare the original ++ event rates to — — rates, 
with both beams flipped, which should be same as ++ rates if parity is conserved. 

If the luminosity is high enough such that there is a non-negligible probability of 
multiple collisions per bunch crossing occurring, it may present a problem for both 
double- and single-spin asymmetry measurements. If the occurrence of multiple 
collisions per crossing is not spin- dependent, the relative luminosity measurement is 
not affected. If the probability of multiple collisions per crossing is spin-dependent 
due to spin-dependent beam conditions, it can likely be studied using the spin flipper. 
If it is spin- dependent due to a physics asymmetry, the problem is similar to the 
one discussed above in the case of no more than a single interaction per crossing. 

It should be noted that the analogous formula to Eq. 15. II for longitudinal double- 
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spin asymmetries is Eq. IA.1I 

1 iV ++ - KN + - 
AlL = P X P 2 7V++ + KN+- {AA) 
where Pi and P 2 are the polarizations of the two beams, iV ++ represents the par- 
ticle yield from same-helicity crossings (++ and ), N + ~ represents the particle 

yield from opposite-helicity crossings (H — and — h), and 1Z is the relative lumi- 
nosity between same- and opposite-helicity crossings. For longitudinal double-spin 
asymmetries, because there is no formula analogous to Eq. 15.81 in which luminosity 
differences cancel to several orders, it is essential to determine the relative luminosity 
accurately. 

The transverse double-spin asymmetry, Att, is expected to have the form A + 
Bcos2ip, where A and B are constants. Therefore there is a potentially non-zero 
constant term as well as an azimuthal dependence. The same techniques for un- 
derstanding the relative luminosity as in the double-longitudinal case should be 
applicable. 

In the case of a transverse single-spin asymmetry, A^, as in the present analysis, 
the form is expected to be purely azimuthal, expressible as A cos (p. In this case 
there is no concern regarding physics asymmetries present in the detectors used as 
the MB trigger and to measure the relative luminosity, as long as the detectors cover 
the full 2ir in azimuth. A new concern arises, however. An azimuthal dependence 
of the detector efficiency could produce bias in the "minimum-bias" trigger and the 
relative luminosity measurement. Comparing the asymmetry results obtained using 
Eq. 15.81 which is largely insensitive to the relative luminosity, provides a handle on 
this issue. In addition, many of the cross checks and systematic studies discussed 
above for the double-spin asymmetries would also be relevant. 



